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Abstract. We restate Kay's 1998 hypothesis which simultaneously offers an objective 
definition for the entropy of a closed system, a microscopic foundation for the Second 
Law, a resolution of the Information Loss (and other) Black-Hole Puzzle(s) as well as 
an objective mechanism for decoherence. The hypothesis presupposes a conventional 
unitary theory of low-energy quantum gravity and offers all this by taking the physical 
density operator of a closed system to be the partial trace of its total density 
operator (assumed pure) over gravity and by defining its physical entropy to be its 
'matter-gravity entanglement entropy'. We further recall Kay's 1998 modification of 
non-relativistic (many-body) quantum mechanics based on Kay's hypothesis with a 
Newtonian approximation to quantum gravity. In this modification, we find formal 
expectation values for certain 'observables' such as momentum-squared and Parity 
are altered but those for functions of positions are unaltered. However, by arguing 
that every real measurement can ultimately be taken to be a position measurement, 
we prove that, in practice, it is impossible to detect any alteration at all and, in 
particular, we predict no alteration for an experiment recently proposed by Roger 
Penrose. Nevertheless, Kay's modification contains no Schrodinger Cat-like states, and 
also allows a possible 'events' interpretation which we tentatively propose and begin 
to explore. We also obtain a first result of Second-Law type for a non-relativistic toy- 
model closed system and argue that similar results will apply for a wide class of model 
Newtonian and post-Newtonian closed systems although we argue that ordinary actual 
laboratory-sized systems can never be treated as closed for the purpose of calculating 
their entropy. Compared with 'collapse models' such as GRW, Kay's Newtonian theory 
does a similar job while being free from ad hoc assumptions. 
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1. Introduction 

1.1. Background 

In 1998, [1], one of us (BSK) argued that several of the problems and puzzles 
concerning decoherence and thermodynamic behaviour, which are inherent in our 
current understanding of quantum physics in general, as well as several problems and 
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puzzles specifically related to quantum black holes would appear to find natural 

resolutions if one makes the following hypothesis (which is in three parts): First, that 
for the resolution of these problems and puzzles, a quantum gravitational setting is 
required; more precisely, a low-energy quantum gravitational setting, where, 'low' here 
refers to energies well below the Planck energy. Second, that low energy quantum gravity 
is a quantum theory of a conventional type with a given closed quantum gravitational 
system described by a total Hilbert space H tota j which arises as a tensor product 

^total = "^matter ® ^gravity (1) 

of a matter and a gravity Hilbert space <EI1>, and with a total time-evolution which is 
unitary, but that, while, as would usually be assumed in a standard quantum theory, 
one still assumes that there is an ever-pure time-evolving 'underlying' state, modelled 
by a density operator of form 

Ptotal = I (2) 

\& G ^totab a ^ eacn 'instant of time', one should add the new assumption that the 
physically relevant density operator is not this underlying density operator, but rather 
its partial trace, p ma tter; over ^gravity Third, that the physical entropy of a closed 
quantum gravitational system is to be identified with the von Neumann entropy of 

Pmatter > 

•^physical = — ^ t r (Pmatter m Pmatter)- 
In other words, <Hvl>: 

The physical entropy of a closed system is its matter- gravity entanglement entropy. 

With this hypothesis, an initial underlying state of a closed quantum gravitational 
system with a low degree of matter-gravity entanglement would be expected to become 
more and more entangled as time increases and thus the physical entropy, as we have 
defined it above, would be expected to increase monotonically, thus (when the theory 
is applied to a model for the universe as a whole) offering the possibility of an objective 
microscopic explanation for the Second Law of Thermodynamics and (when the theory 
is applied to a model closed system consisting of a black hole sitting in an otherwise 
empty universe) offering a resolution to the Information Loss Puzzle. It offers both of 
these things in that, by defining the physical entropy of Ptotal t° be the von-Neumann 
entropy of p ma tter; on e reconciles an underlying unitary time-evolution on T^total w ith 
an entropy which varies/increases in time <TyJ>. 

A second paper, [2] (also by BSK) also in 1998, investigated the implications of 
the hypothesis of [1] for the decoherence of ordinary matter in the 'Newtonian' - i.e. 
non-relativistic, weak-gravitational-field - regime. With some further assumptions, the 
conclusions of this work were that, if, in ordinary non-relativistic quantum mechanics, 
the centre-of-mass degree of freedom of a uniform-mass-density ball of mass M and 

| The small Roman numerals in angle- brackets refer to the section entitled 'Notes' (Section fTOj) at the 
end of this article. 
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radius R would be described by a Schrodinger wave function ip e L 2 (M 3 ) then the 
physically relevant density operator p ma tter (from now on, where it can cause no 
confusion, we will often just call this p) is given by 

p(x,x')=p (x,x')e- D ( x < x '\ (3) 

where p (x, x') = ^(x)?/'*(x / ) is the position-space density matrix one would expect to 
have were gravitationally- induced- decoherence effects to be ignored, and D, which is 
called the decoherence exponent, is a certain function <lvil> of x and x' which depends 
on x and x' only through the square of their difference (x — x') 2 and which vanishes 
when x = x'. Explicitly, D is given (below and from now on, a stands for |x — x'| and 
we regard D as a function of a) by 

D(a) = 216M 2 r W")-^«*(/Q) 2 / na/R - sin( K g/R)\ ^ 
Jo K V Ka / ^ / 

and we remark that this is a monotonically increasing function of a [2]. 

To summarize, the overall effect of the new theory is to replace the 'would- 
be' Schrodinger wave function ip (or more precisely the would-be density operator 
Po = IV'XV'I) by the physical density operator p of (J3j). 

Moreover, in [2], an asymptotic limit of D was obtained for the case where a ^ R, 

D(a) ~ aa 2 + 0(aV \n(a/R)), (4) 

where (by default we shall assume Planck units where G = c = h= k = l)a = 9M 2 /R 2 . 
We refer to P| as the Gaussian asymptotic regime. 

We remark that we expect this Gaussian regime to always give a good 
approximation for a ball with a mass around or bigger than the Planck mass, since 
then, whenever a fails to be very much less than R, e~ D ^ will anyway be very small. 

We also remark that, as long as the ball in question is far from being on the verge 
of collapsing to form a black hole (which we anyway expect it to need to be for our 
Newtonian approximation to be valid) then the constant a in the Gaussian regime (pD) 
will (in Planck units) be very much less than one. 

An asymptotic limit was also obtained in the case where a>i?, 

D{a) ~24M 2 lnQ) +0(1). (5) 

in view of which we refer to this limit as the logarithmic asymptotic regime. 

We shall sometimes study a model one- dimensional analogue of ([3]) which is taken 
to be exactly Gaussian, so that the relation between the physical and would-be density 
operator is given by 

p(x,x') = p (x,x')e-^-^ 2 , (6) 

where 

po(x,x') = ip(x)ip*{x). 

where if) is taken to be a 'would-be' wavefunction on the line and ft is a positive constant. 
We shall call this the one- dimensional Gaussian model. 
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In view of (J3J) and the subsequent remark, we would expect that, if we identify k 
in with a = 9M 2 /R 2 , then, for M around or bigger than the Planck mass, this one- 
dimensional Gaussian model will give a good approximate description of the relation 
between the would-be wavefunction and the physical density operator describing the 
centre-of-mass degree of freedom of a uniform mass bead (i.e. a uniform mass ball with 
a narrow straight tube bored through its middle) of mass M and radius R constrained 
to slide on a straight wire, where x is taken to denote the position of the centre-of-mass 
of the bead along the wire relative to some choice of origin. 

We note that if such a bead is constrained to move between hard stops, say a 
distance 2(R + 5) apart - so that, calling the position of the centre of the bead when 
it is half-way between the stops 'the origin', any would-be wave function, ip(x) for the 
centre-of-mass motion is required to vanish at x — ±<5 - then the Gaussian model would, 
again in view of (jlj), be expected to give a good approximation for any mass M provided 
5<R. 

[2] also generalized the formula ([3]) to the many-body wave function (for the centre 
of masses) of a system of many balls. The general formula, for a system of N balls takes 
the form 

p( Xl , . . . , Xjv ; x' 1; . . . , x'jv) = p (xi, • • • , x^; x^, . . . , X ' jv ) e -D(xi,.,x^;xi,...,x^ v ) ) ( 7 ) 

where p (xi, ..,xn;x[, ...,x' N ) takes the form \I/(xi, .., xjy)^*(x[, x' N ) for some many- 
real- variable wave function \1/ and D is a certain function which depends on xi . . . xjv 
and ~x! x . . . x'jy only through the distances |xj — Xj|, |x' 7 — xj|, |x/ — x'j| and |x' 7 — x'j| 
(J, J = 1 . . . N) and which vanishes whenever xj = x 7 for all I = 1 ... N. 

In much of this paper, we shall focus our attention on states involving a single ball 
in 3 dimensional space or on our even simpler one-dimensional Gaussian model but we 
note that in the generalization to the case of N balls of equal masses and radii, there 
is a generalization of the Gaussian asymptotic regime in the leading term of which D 
takes the explicit form <dviib> 



D(xi, . . . , Xjv; x' 1; . . . , x.' N ) = a (xi + x 2 + ... + x w - x[ - x' 2 - ... - x 



N i 



This leads us to define an iV-body generalization of the one-dimensional Gaussian model 
by 

P ( Xl , .., X N] X[, X' N ) = po{ Xl , .., X N] X[, x ' jv ) e -«(xi+x 2+ ...+x„-x 1 '-x 2 '-...-x ,v') 2 

where again po(xi, .., xjv; x[, x' N ) takes the form *$?(xi, .., xat)\1/*(x / 1 , x' N ) for some 
many-real- variable would-be wave function We shall refer to this latter for one 
of the results obtained in the Appendix (Section [9]). (Note <fvlib> also includes the 
generalization to the many-ball case of the logarithmic asymptotic regime.) 

[2] further proposed a Schrodinger picture time-evolution rule appropriate to this 
Newtonian approximation, according to which (in the single-ball case) the time-evolving 
physically relevant density operator is obtained by adopting the formula ([3]) at each 
moment in time, t, so that 

p(t;x,x')=p (t;x,x')e- D ( x ' x '), (8) 
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where po(t; x, x') = ^(x)^* (x') where the would-be wave function ip t evolves exactly as 
in ordinary non-relativistic quantum mechanics, for some choice of Hamiltonian, H, so 
that V ; t( x ) = (x\e~' lHt \ijjo) < lviiil >- (Similarly, in the many-ball case and the Gaussian 
model etc., we adopt formula (JTj), respectively (jUJ) etc. at each moment in time, again 
assuming that the relevant would-be wave functions evolve in time according to some 
choice of Hamiltonian which we again call H.) 

For a Schrodinger Hamiltonian of form H = —V 2 /2M + V, it may easily be checked 
that this amounts to the statement that p(t; x, x') evolves according to the master 
equation 

P = ^ (V 2 - V' 2 + V(x) - VM) p-^V a D{x-x') (-iV a - iV a )p.(9) 

We remark that it is easy to see that the analogous master equation for the one- 
dimensional Gaussian model can be written in the simple operator form 

p=-i[H,p]-^[x,\p,p]]. (10) 

While these master equations can be shown not to have 'Gorini-Kossakowski- 
Sudarshan/Lindblad' (cf. [6], [7]) form (we shall abbreviate this elsewhere to 'GKS/L 
form') and even not to arise from a time-indexed family (i.e. not necessarily semigroup) 
of 'completely positive maps' acting on an initial p (see <Hxb> for an explanation of all 
this) they clearly do have solutions, specified by the formulae (JED, (EJ), which consist of 
a density operator <|xj> p t at each instant of time and (in either case) the map from 
Pt to p t+s arises in the form p t+s = p{s)pt where the operators (on the space of density 
operators on our matter Hilbert space) p(s) clearly form a one-parameter group (not 
just semi-group). 

1.2. Purpose of present paper and summary of results 

The project which gave rise to the present paper had three purposes: Firstly to 
explore the prospects for experimentally testing the theory of PQ, secondly to explore 
whether /how the theory may form a basis for a solution to some of the conceptual 
problems of quantum mechanics (i.e. of the 'measurement problem'), and, thirdly and 
lastly, to explore in more detail what the theory predicts about the time- dependence of 
entropy. 

In the present paper, as a first step, we shall begin to explore all three of these 
issues, proceeding as if, for the relevant systems, the Newtonian approximation of [2] 
applies. It was argued in [2] that this assumption is unjustified for physically realistic 
laboratory-sized systems, at least as far as questions of entropy are concerned - the 
problem being that, while such systems will satisfy the necessary approximations in the 
sense that the relevant speeds are slow enough and the relevant mass distributions very 
far away from black-hole collapse etc., it is not expected (for the reasons listed in [2] 
and discussed further in Section [Zj) to be appropriate to treat them as closed. (They 
should rather be treated as open systems as indicated in Note <EiH>.) 
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In this article, the systems we treat will nevertheless be treated as if they were 
closed. Thus considerable caution should be exercised in assigning a direct physical 
interpretation to our results. However, we shall argue in the later parts of our Discussion 
section (Section [7]) that, as far as our first two questions (experiments and measurement 
problem) are concerned, the answers we obtain will, when properly interpreted as we 
explain in Section [H survive when we drop this assumption and that, as far as our third 
question (entropy) is concerned, the results we obtain will at least survive as a part - 
albeit perhaps a small part - of the full results for actual (presumably partly relativistic) 
closed systems and also be of interest as providing 'toy models' for these full results. 

Proceeding then as if the relevant systems were closed and our Newtonian theory 
of [2] were exact, we amplify on each of these issues in turn: [2] pointed out that, as 
far as the mathematical formulae are concerned, the Newtonian theory of [2] leads to 
large effects in the 'amount of decoherence' already for 'small macroscopic' systems. 
Thus e.g. a ball with 10 times the Planck mass with a would-be wavefunction for its 
centre-of-mass a Schrodinger Cat-like superposition 

tp = + C 2 1p 2 (11) 

(cf. (1671) in Note <lvil>) where ipi and ip2 are sharply localized around two distinct 
locations separated by a distance a equal to one tenth of the ball's radius, will, by 
01]), have a decoherence exponent, -D(a), of approximately 9, and thus have a physical 
density operator very close to that for a statistical mixture 

p=| Cl |Vl)(^l| + |c 2 |V2)^2| (12) 

of states localized at the two locations. More precisely p will be as in fl68|) with {g\\g2} 
approximately equal to e~ 9 . It thus may seem at first sight plausible that, with suitable 
experimental setups, differences from the predictions of standard quantum mechanics 
should be experimentally detectable in the laboratory. Moreover, there have been 
several proposals (the proposals we mention below are [8], [9]) for such experimental 
setups, albeit the theories which the authors of these setups had in mind were somewhat 
different from the theory of [I], [2]. In order to explore what the prospects are for such 
experimental tests, we need of course to take a view about how our mathematical 
formalism relates to experiment. To this end, we have explored the consequences 
of adopting what we shall call here the naive pragmatic interpretation according to 
which the usual sorts of 'observables', e.g. (for the single ball model) position (x) and 
momentum (p) and functions of these, parity (V) and so on, are represented by the 
usual self-adjoint operators and the expectation value of the observable represented by 
the self-adjoint operator A obtained in the state p is given by the usual formula 

The expectation value, (A) p , of A in the state p is equal to tr(pA). (13) 

As we report in Section [2] we find that, with this naive pragmatic interpretation, the 
expectation values of x and p and also of arbitrary functions of x are unchanged, 

tr(px) = tr(p x), tr(pp) = tr(p p), tr(p/(x)) = tr(p /(x)), (14) 
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but that expectation values for p 2 and parity, V, acquire extra terms (see (1161) and 
(j!8p ). In Section [31 we explore the change in the naive expectation values for p 2 further 
and point out that these lead to a non-zero state-independent minimum value for the 
expectation value of p 2 and, related to this, to a modification of the usual x p Heisenberg 
uncertainty relation. (Section [3] stands to one side of the main line of the paper. The 
results are expected to be of interest in their own right but are not needed or referred 
to in subsequent sections.) 

We shall however (i.e. in spite of the alterations in naive expectation values of 
some observables from the those of standard quantum mechanics) argue that our naive 
pragmatic interpretation is too naive. As we shall discuss in detail in Section HI if one 
considers an experimental setup e.g. to determine the deviation of the expectation value 
of p 2 for a ball in a particular physical state from its expectation value according to 
standard quantum mechanics in the corresponding would-be wave function, one would 
presumably need to introduce at least a second 'probe' particle and a measurement of 
p 2 for the original ball might then e.g. amount to some sort of position measurement 
involving the probe particle. However, as long as our probe particle is also non- 
relativistic <dxib>. the two-body versions of our formulae (TT3I) are assumed to apply 
and one can easily see from these that they entail, when formula (fT3l) is applied at the 
two-body level, that 

Any sort of position measurement on the relevant physical density operator will always 
give the same result as that predicted, according to standard quantum mechanics, for the 
same measurement on the corresponding would-be wave function. 

In the case discussed above, the relevant physical density operator is the two-body 
density operator replacing the time-evolute of the tensor product of the would-be wave 
function for the ball and the initial would-be wave function for the probe particle. We 
will prove a generalization of this result for iV-bodies in Section H] and we shall refer to 
this latter generalization as our Position Measurement Theorem. 

One might of course think of some cleverer setup in which the measurement on the 
probe particle is also not a sort of position measurement, but we shall advocate the view, 
that for a measurement to have been convincingly made, then, sufficiently far along the 
Heisenberg- von Neumann chain of measurements (see Chapter 6 in [TO]) it eventually 
ends up, or can be regarded as ending up, as a position measurement - traditionally 
one talks of the position of a 'pointer' on a 'dial'. We shall call this view the corrected 
pragmatic interpretation; it is still 'pragmatic' insofar as it still applies the formula ffl3l) 
at some stage along the Heisenberg- von Neumann chain. Adopting this interpretation, 
the general A-body argument we mentioned above applies and we thus have the rather 
surprising conclusion: 

On our corrected pragmatic interpretation, the experimental predictions of the 
Newtonian approximate theory of JE/ are identical with those of standard quantum 
mechanics! 

(A) 
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Of course we expect that, in the fully relativistic strong- field regime, the hypothesis 
of [lj will predict (large) differences but we expect them to depend on the post- 
Newtonian aspects of the theory. In Section H] (see especially Note < lxxvil >) we give 
an estimate which suggests that these post-Newtonian corrections will, however, be 
completely negligible for the sort of laboratory experiments which are currently under 
consideration. 

Of particular interest in connection with our conclusion (A) above is the photon- 
interferometer-based experiment [21], [22] proposed by Roger Penrose to detect 
modifications to the predictions of standard quantum mechanics predicted by theories 
of 'collapse-model' type such as that of Ghirardi- Rimini- Weber [13] (elsewhere in the 
paper, we refer to this as the 'GRW model) and others and, in particular, by the 
gravity-induced quantum state reduction proposal [20] [21] [22] of Penrose himself. At 
the heart of the experiment is a macroscopic object (crystal) attached to a movable 
interferometer mirror which, by the arrangement of the experiment would get put, 
according to a standard quantum-mechanics analysis, into a superposition, one branch of 
which remains at rest, while the other oscillates between two distinct locations. If (e.g. 
because of Penrose's gravity-induced quantum state reduction) the crystal decoheres 
as a result of visiting (i.e. in one of the branches of the superposition) the distinct 
location during its oscillation cycle, a superposition of the photon state in the two 
arms of the interferometer will (by the design of the experiment - see Section 14. If) get 
converted into a partly decohered mixture thereby changing the predicted exit-paths of 
a certain fraction of the photons travelling through the interferometer from what would 
be predicted by standard quantum mechanics, and this change can then be detected 
with suitable photon detectors. 

According to Kay's theory [2] such a crystal will (if the size is sufficient and 
the distinct location suitably distant) in some ways similarly to what is predicted by 
Penrose's quantum state reduction, be predicted to decohere each time it visits its 
distinct location and yet, according to our conclusion (A) (and assuming that it applies 
when some of the particles involved are photons <lxil>) what is detected at the detector 
must be identical to what would be predicted by standard quantum mechanics. 

In Subsection 14.11 we explain how this difference in predictions between the theory 
of Kay and the quantum state reduction of Penrose may be traced to the following 
qualitative difference between their proposed decoherence mechanisms. On the Penrose 
proposal it is (at least implicitly) assumed that once the crystal has decohered, then it 
cannot later recohere. In the Kay theory of [2J, on the other hand, the superposition 
involving the crystal is understood to decohere when, in one of the branches of its 
superposition, the crystal is in its distinct location, but to recohere each time the crystal 
returns to being in the same location in both branches. 

Turning to the second of our purposes, which, we recall, was to explore how the 
theory of [1] may form a basis for a solution to some of the conceptual problems 
of quantum mechanics, we begin by remarking that while the corrected pragmatic 
interpretation described above seems the best we can do at the moment to interpret 
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our theory and make experimental predictions, it is still unsatisfactory in that, like the 
standard interpretation of standard quantum mechanics, it makes essential reference to 
'observables' and hence implicitly to 'observers' and, as e.g. John Bell eloquently argued 
(see e.g. [12]) one might hope that our current understanding of quantum mechanics 
will one day be superceded by a theory which is more objective in nature. Indeed 
many authors have asked for a modification of quantum mechanics (or maybe for a 
changed interpretation of standard quantum mechanics) in which there is an objective 
notion of 'events' which 'happen'. For theories which are couched in terms of a time- 
evolving density operator p(t), as in either the full hypothesis of [1] or the Newtonian 
approximate theory of [2] being considered here (but also 'collapse models' which we 
mentioned above in connection with the Penrose experiment)) a natural proposal would 
seem to be < jxivb > <3xvl> < Jxvil > : 

The set of possible events which can happen at a given time, t, is to be identified with 
the set of spectral subspaces of pit) at that time. The probability with which a given such 
event occurs is to be identified with m\ where A is the eigenvalue of pit) belonging to 
that subspace and m is its multiplicity (i.e. the dimension of the subspace). 

(B) 

As we will discuss further in Section El this proposed interpretation would, for 
example, interpret the state (JT2"1) in terms of two possible events, namely the subspace 
spanned by ipi and the subspace spanned by tp 2 , and assign to these events the 
probabilities |ci| 2 , | C2 1 2 respectively. We shall also explore, in Section El what this 
proposed interpretation entails for the Gaussian model (jSJ) and shall actually simplify 
this further by studying that only to first order in k. 

In order to do this, we obtain, in the Appendix, the diagonalization of the density 
operator of the one-dimensional Gaussian model to first order in k. A typical result is 
that, for a given would-be wave function ip which is an even function or an odd function, 
the density operator p in ([6]) takes, to first order in k, the form 

p = Ai|0i)(0i| + A 2 |0 2 )(02| 

where Ai = 1 — 2k{i/j\x 2 iP) is (for small k) close to 1 and <\>\ (see (|33|) ) resembles ip while 
A 2 = 2k((iI)\x 2 iP) is small and 02 (equal to a normalization constant times xip) has the 
opposite parity to ip. According to our proposal, this would be interpreted to mean 
that, for such an (even or odd) would-be symmetric quantum wavefunction on the real 
line, there are two events, one with probability Ai which is identified with the subspace 
spanned by 0i and one with probability A 2 which is identified with the subspace spanned 
by 4>2- 

Concerning our third purpose, which is to explore what the theory of [1] has to say 
about the time-evolution of entropy, we make a small beginning in Section [6] by exploring 
how the von- Neumann entropy S and also the entropy-like quantity Si = 1 — trp 2 
vary in time, mainly treating our one-dimensional Gaussian model restricted further 
to would-be wave functions which are either even or odd and to the regime where 
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working to first order in k suffices so that we may utilize the diagonalization results 
obtained in the Appendix. As is apparent from the formulae quoted above for the 
eigenvalues Ai and A2, the values of S and S\ for this model when the system is in 
a given would-be state, if), at a given time t, will be correlated with the value of the 
standard-quantum-mechanical quantity { r ip{t)\x 2 'ip{t)) at that time, which, in standard 
quantum mechanics, is of course a measure of the 'spreading of the (would-be) wave- 
packet'. This leads, if, as is the case for a free Hamiltonian, the Hamiltonian governing 
the would-be dynamics leads to spreading of the wave packet at times on either side of 
a time of minimum spread, to a 'two-sided entropy increase' result. On the other hand, 
for a would-be Hamiltonian, say, of harmonic-oscillator type, the would-be wavefunction 
will be periodic in time, and therefore the amount of wavepacket-spreading and therefore 
also the values of S and S\ will be periodic in time, indicating a continual process of 
decoherence and recoherence. We point out that this is analogous to the decoherence- 
followed-by-recoherence mentioned above in our discussion of the Penrose experiment. 
We also compare and contrast this Hamiltonian-dependent behaviour of S and S\ 
with the typical sort of monotonic entropy-increase result that one finds for density 
operators which satisfy master equations of GKS /L form which are usually the relevant 
master equations in collapse models. Specifically, we compare and contrast the results 
mentioned above which are relevant to the master equation ({TO]) with an easily- derived 
result on the Hamiltonian-independent monotonic increase of entropy for the Barchielli- 
Lanz-Prosperi [27] master equation fl38l) 

P = -i[H,p) - c[x, [x,p]) 

(c a positive constant) which may be regarded as a prototype master equation of GKS/L 
form. (Elsewhere in this paper, we shall refer to this as the 'BLP' equation.) 

The main part of the paper ends, in section [7J with a discussion of various issues 
arising from the earlier sections. 

We shall summarize here some of the main overall conclusions which we draw, in 
Section [7J by comparing and contrasting what we have learned about the theory of [2] 
with the properties of 'collapse models' such as GRW [13]. First, we conclude that, 
for ordinary laboratory-sized quantum systems, the cautionary remarks mentioned at 
the start of this subsection can be ignored as far as quantum mechanical measurements 
are concerned and (in view of our Position Measurement Theorem and our 'corrected 
pragmatic interpretation') our theory really does have identical experimental predictions 
to standard quantum mechanics (although there will be deviations when relativistic 
effects are taken into account) while, at the same time, being free from macroscopic 
Schrodinger Cat-like superpositions. This is to be compared and contrasted with collapse 
models which predict small deviations from the predictions of standard (non-relativistic) 
quantum mechanics while also being free from Schrodinger Cat-like superpositions. 
Another feature of collapse models is that they give rise to an increasing amount 
of decoherence as time increases. For our theory, we give, in Section [7J plausibility 
arguments for a two-sided entropy- increase result, of the sort discussed in Section 
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E] and mentioned above, for a wide variety of 'generic' model closed systems within 
the Newtonian framework of [2] and we also argue for a plausible post-Newtonian 
extension of such results. Concerning the two-sidedness, we advocate, in Section 
[7J the natural physical interpretation that the negative-times should be discarded as 
physically irrelevant and the time-zero density operator regarded as the 'initial state'. 
(If the closed system is taken as modelling the universe, then this will model the initial 
physically relevant density operator of the universe.) In this sense, such results would 
amount to proofs, for appropriate models, of the Second Law. As far as decoherence 
is concerned, such Second-Law results can be interpreted as entailing an increasing 
amount of decoherence as time increases. However, we do not interpret these results 
as applying to actual ordinary (non-relativistic) systems of laboratory size because 
indeed we will argue in Section [7] that, as far as questions of entropy are concerned, 
actual closed systems must be partly relativistic and indeed may well be of galaxy- 
size! Rather, we argue that actual small ordinary laboratory-sized systems need to be 
regarded as open systems (See Note <lxiil>) and that the correct mechanism behind their 
decoherence is essentially the same as that advocated on the traditional environment- 
induced decoherence paradigm. (The entropy, i.e. amount of decoherence, will then 
increase or decrease with time according to the model and state etc.) (We have omitted 
here a comparison of approaches to 'events' in the two theories. For this, we refer to 
Sections [5] and [7] and especially Note <Jxvl>.) 

Our overall conclusion from this comparison is that our approximate Newtonian 
theory of [2] does a similar job to the job which 'collapse models' such as GRW [32] do 
- albeit in a slightly different and partly unexpected and surprising way but with the 
advantage that, while collapse models are ad hoc, our Newtonian approximate theory 
is part and parcel of a general hypothesis (i.e. the hypothesis of [T]) which also resolves 
a number of other puzzles. 

This paper arose out of an ongoing research project of one of the authors (BSK). 
While both authors were involved with the new results presented in the paper, the main 
involvement of VA was with Sections EJ [3], and M with Section [3] and Subsections 19.11 and 
I9.2l being mainly due to VA. Some further material related to these sections can be found 
in [4]. Some parts of the rest of the paper (including the last paragraph of Subsection 
11.21 the proposed interpretation in terms of events in Sections [1] and [5] and the Position 
Measurement Theorem discussed in Sections [1] and HI including the discussion of the 
Penrose experiment in Subsection 14. ip ) are mainly due to BSK, as are Sections [6] and [3 

The paper concludes with an extensive Notes section (Section [10l due to BSK) 
which was needed to put the results of the main body of the paper in context and 
clarify their significance. It both collects together some relevant background material 
and also includes some new material which is expected to be of interest in its own right. 

In particular, aside from their importance as relevant background to the new results 
reported in this paper, the first paragraph of Subsection 11.11 together with Notes <0>, 
<dmt>. <Hvt>. <]v> and <jxiib> will, we hope, be of interest as constituting a self- 
contained, updated (i.e. in the light of more recent work of BSK) re-statement of 
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the hypothesis, which was first made in [T], for a resolution of the puzzles listed in Note 
<S>. 

2. Naive pragmatic interpretation: trace formulae 

It is easy to see from ([3]) that, if A is any function, /(x), of position and also if A is the 
momentum operator p, then tr(pA) is equal to tr(p A). To see the former, notice that 



tr(p/(x)) = J /(x)p(x,x)d 3 x 

and that, by virtue of Z?(x, x) = 0, 

p(x,x) = p (x,x). 

For the latter, we combine 

f d , 

fc(PPa)= / -«^P(x,x')|x=x'(i 3 X / 



and 

._d , . ( 

dx 



p(x,x') = -i (A po ( XjX ^ e-DW) _ 2mpo (x,x')(x a 
+higher order terms which vanish when x = x' 



(here we used (J1J) and notice that all but the first term in the second equation will 
vanish when x = x'. 

On the other hand, one easily sees ('a' here refers to a single index and is not 
summed over) from 

-^p(x,x') = - (^Po(x,x')) e-^> - 4a (^Po(x,x')) (*" - x'^ D ^ 

+2apoe~ D( ^ x ' x - ) + higher order terms which vanish when x = x' 

that, for each single Cartesian component of momentum squared, 

ti(pp 2 a ) = tr(p p2) + 2a, (15) 

whereupon, summing over a, we clearly have 

tr(pp 2 ) = tr(p p 2 ) + 6a. (16) 

It is interesting to note that all the above results remain unchanged if one replaces 
the term e _I? ( x ' x ') in (jHJ) by its Gaussian asymptotic form e ~ a ( x ~ x ') 2 or even if one replaces 
the latter by its 'first-order in a 1 form 1 — a(x — x') 2 , albeit the resulting expression 
for p is then not to be taken seriously except when x and x' are sufficiently nearby for 
a(x — x') 2 to be small compared to 1. In this spirit, we write 

Papprox(x, x) = p (x, x) (l - a(x - x') 2 ) 

i.e. 

Papprox = po ~ ax a x a p + 2ax a p x a - ap x a x a = p - a[x a , [x a , p ]] 
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(summed) whereupon we have 

tr(papprox^) = tr(p Q A) - atr([x a , [x a , Po ]]A) = tr( Po A) - atr{p [x a , [x a , A}]) (17) 

(summed) where we have used the easy standard trace identity, tr(A[B,C]) = 
tr([A, B]C), etc. in the last equality. One easily checks that this latter equation 
reproduces the exact expectation values for /(x), p and p 2 (unsummed) (and p 2 ) found 
above. 

On the 'naive pragmatic interpretation' discussed in the introduction, the equation 
tr(p/(x)) = tr(p /(x)) predicts e.g. an identical diffraction pattern to that of standard 
quantum mechanics e.g. for a double-slit experiment involving a ball of ordinary matter 
even as heavy as or heavier than the Planck mass (cf. the experiments of Zeilinger et al 
[T9] albeit for smaller masses). On the other hand, if we consider the wavefunction for, 
say, the ground state of such a ball in a spherical or cubical box etc., then the results 
( |T5l) . (TIB]) would seem to predict a detectable difference in the expectation value of the 
squared momentum. 

One can also calculate tr(pV) where V is the parity operator. 

tv(pV) = y"p (-x,x)e- D( - x ' x) rf 3 x 

Assuming the above 'first-order in a' approximation to the Gaussian approximation 
held (this would be the case e.g. for the centre-of-mass wavefunction of a ball with mass 
much bigger than the Planck mass and radius R confined to a cubical or spherical box 
of diameter not much bigger than 2R) we would have, to first order in a 

tr(pV) = tr(> P) - atr(p [x a , [x\V]}) = tr(p P) - 4a(^|x 2 P^)- (18) 

We remark that, if the would-be wave function ?/>(x) has even parity, then this is 
1 — 4a(?/>|x 2 '?/>) i.e. 1 minus 4a times the squared uncertainty in x of the would-be 
wave function and, if t/>(x) nas °dd parity, then it is —1 + 4a(?/>|x 2 ?/>) i.e. —1 plus 4a 
times the squared uncertainty in x of the would-be wave function. Also, the counterpart 
to this result for our one- dimensional Gaussian model (J6j), is that, for a would-be wave 
function ip{x) on the line 

tr(pV) = tr(p V) - Ktr(p [x, [x, V}}) = tr(p V) - 4/t(^|x 2 P^). (19) 

This would be appropriate e.g. to the bead-on-a-wire constrained to move between 
sufficiently close-together hard stops as discussed in Subsection 11.11 Moreover, if the 
wave function is even, we'll have 

tr(pP) = 1 -4a(^|x 2 ^) (20) 

(and if it is odd, 

tr(pP) = -1 + 4a(^|xV)). (21) 

We will discuss the significance of ([TBI . (JT9~]1 . ff20l and f|2TT) further in Section [51 
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3. Minimum Ap 2 and modified uncertainty relations 

(Note The results of this section are not referred to or needed in the remainder of the 
paper. Note also that, throughout this section, we shall not assume that a repeated 
index is summed over.) 

Continuing to adopt the naive pragmatic interpretation, it is interesting to note 
that Eq. ([15]) implies that there is a minimum value (i.e. 2a = 18M 2 /R 2 ) for the 
expectation value (p 2 ) = tr(pp 2 ) of the square of each component of the momentum in 
a given physical state for the centre of mass motion of one of our balls of mass M and 
radius R. Given that tr(pp a ) = tr(p p a ), and defining the squared uncertainty, (Ap a ) 2 , 
in p a in the usual way by 

(Ap a ) 2 = tr(pp 2 a ) - tr(pp a ) 2 , 

we have 

(Ap a ) 2 = tr(p p 2 ) - tr(p p a ) 2 + 2a = (A oPa ) 2 + 2a (22) 

where A p a is the usual uncertainty in p a of the corresponding would-be wave function. 
So the squared uncertainty in each component of momentum also has a minimum 
possible value of 2a. 

To get some idea of the orders of magnitude involved, we can convert the minimum 
expectation value, 6a, of p 2 into a 'maximum wavelength' (cf. the de Broglie relation 
'A = 27i7i/|p|'), so that 




2nR 
3V6M' 



In any system of units, this is (27r/3v^6)(Mp^ anc ] t /M)_R. Thus, for the centre-of-mass 
motion of a uniform density ball of mass around the Planck mass (i.e. around 10~ 5 g) 
A will be of the same order of magnitude as the ball's radius. For a proton, modelled as 
a uniform density ball with the proton mass and a radius around 10~ 13 cm, A will be of 
the order of 10 kilometres! 

In view of tr(p/(x)) = tr(p /(x)) (and with obvious definitions) the squared 
uncertainty (Axb) 2 in a given physical state for a given component of position will equal 
the usual squared uncertainty (A Xb) 2 in the corresponding would-be wave function. 
Combining this with (1221) leads to the replacement of the usual Heisenberg uncertainty 
relation 

(A x a ) 2 (A p fe ) 2 > -S ab 

by 

(Ax a ) 2 (Ap fe ) 2 = (A x a ) 2 ((A p fe ) 2 + 2a) > ~6 ab + 2a {Ax a f . (23) 

We remark that if we were to take this relation to be fundamental then we could 
recover the lower bound Ap a > 2a from it. 
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It is interesting to notice that, in the past few years, and partly motivated by string 
theory, many authors have proposed modified uncertainty relations along the schematic 
lines 



where c is a constant. Sometimes these are thought of as applying to spacetime 
coordinates but sometimes to the coordinates of a non-relativistic particle, see e.g. [18J. 
In the latter case, these resemble (12111) - and even more closely, in the case a = b, its 
approximation to first order in a: 



except that the roles of Ax a and Ap a are interchanged, so that, e.g. in [18] and other 
related references, one deduces a lower bound on Ax rather than on Ap. 

4. The experimental indistinguishability from standard quantum mechanics 
and the corrected pragmatic interpretation 

In this section, we elucidate further and prove, the Position Measurement Theorem 
which was outlined in the introduction, according to which any sort of position 
measurement in the Newtonian approximate theory of [2] based on the hypothesis of 
PQ will, under the naive pragmatic interpretation, have an outcome indistinguishable 
from that predicted by standard quantum mechanics. We also amplify on what was said 
in the introduction about the corrected pragmatic interpretation and on the conclusion 
that, with this interpretation, the Newtonian approximation of [2] is experimentally 
indistinguishable from standard quantum mechanics. We also discuss in detail in 
subsection 14.11 how this works out for the class of experiments recently proposed by 
Roger Penrose which we mentioned in the Introduction. 

Consider first our example (introduced in Section [T] where we pointed out that it is 
approximately described - when 8 -C R — by the one-dimensional Gaussian model, also 
introduced there) involving a uniform mass-density bead of mass M and radius R free 
to slide on a straight wire which terminates at stops a distance 2(R + 5) apart. Calling 
the direction of the wire 'the x-axis' and the position of the centre of the bead when 
it is half-way between the stops 'the origin', the would-be wave function, ip(x), of the 
centre-of-mass degree of freedom of the bead in its ground state will, of course, be given 



and thus according to standard quantum mechanics, its momentum in the x-direction 
will take one of the two values ±7r/2<5 with equal probabilities and hence the expectation 
value, (p 2 ), will be 7r 2 /45 2 , whereas, by ( |T5l) . according to our naive pragmatic 
interpretation, it will be approximately 





by 




(24) 




(25) 
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In standard quantum mechanics, one way to experimentally determine the 
expectation value of p 2 (or indeed to obtain the full statistical distribution of p values) 
of the state of such a system might be to rigidly attach a light flat mirror (with diameter 
much larger than S) to the bead, perpendicular to the x direction and, after each time 
resetting the bead to be in the state under investigation, repeatedly to reflect off the 
mirror a small probe particle (with mass much less than M) whose quantum mechanical 
wavepacket is each time prepared to be an approximately monochromatic pencil-shaped 
wavepacket with diameter much bigger than 5, aimed at the mirror along a line at an 
angle of — 7r/4 to the x-axis (see Figure 1), with approximate momentum P pointing in 
the direction of the pencil and say < jxviib > much larger in magnitude than it/ (25). If the 
mirror was fixed rigidly, this pencil would of course turn through a right-angle when it 
reflects off the mirror but because the mirror is fixed to the bead, its (standard quantum- 
mechanical) state will be a superposition (or more generally a mixture of superpositions) 
of wave functions ip n (x) = S^ 1 ^ 2 sm(mrx/25) each of which, in its turn, is of course a 
superposition of states with momenta ±p n = ±mr/2S. Therefore (see Note < jxviil >) 
the probe particle would be predicted to emerge in one of the (approximate) directions 
+7r/4 ± \f2p n jP with a probability depending on the state in the usual way 




Figure 1. Bead on wire with stops fitted with mirror reflecting a pencil of light 

By observing the intensity distribution of probe particles on a suitable screen one 
could thus infer the probability distribution of the state of the bead and, in particular, 
deduce the value of (p 2 )- In particular, if the mirror were in its ground state (T2~4"l) . one 
would have probability 1/2 of the probe particle emerging in each of the (approximate) 
directions +7r/4 ± \p2p\jP (two spots on the screen of equal intensity) and infer 
< p 2 >= p\ = n 2 /{25) 2 while, in any other state, other directions would be populated 
(two more widely spaced spots or more than two spots on the screen etc.) and the value 
inferred for (p 2 ) would be the appropriate standard quantum mechanical value. 

However, we will now argue that, on the Newtonian approximation of [2], according 
to our naive pragmatic interpretation - applied now to a model closed system consisting 
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of our bead-with-mirror plus probe-particle system - the result of such an experiment 
will always be identical with the result predicted by standard quantum mechanics, i.e. 
the result which, according to standard quantum mechanics, would lead one to infer 
that the value of (p 2 ) is the standard quantum mechanical value. Thus, in particular, 
if we always reset the bead to be in its ground state, the value inferred for (p 2 ) will be 
7r 2 /4<5 2 and not the value given in fl25l) . 

To see this, it is convenient to adopt a 'time- dependent' approach to the reflection 
of our probe particles from the mirror. Denoting the ground state of the bead by ip(x) 
and the would-be 'in' wave function for a wave-packet description of the state of one of 
our probe particles at a time well before it hits the mirror by <f)[ n {y), so that the total 
'in' wave function of the bead-probe-particle system is given by 

$ m(^y) = ^(X>0in(y)' 

let us denote by $ Q ut ( x ? Y ) the time-evolute, according to the appropriate standard 
Schrodinger time-evolution, of $j n at the time of observation of the probe particle 
on the screen and denote by p? n the projector |$i n )($i n | and by p° ut the projector 
l^outX^outl- Then the Newtonian approximation of j2] tells us that the total physical 
density operators p- in and p ou t at the two times will be given by (cf. ([7j)) 

p- m (x, y; x, y') = p° n (x, y; x' , y') exp(-D(x, y; x' , y')) (26) 

and 

PoutO^ y; J, y') = Po U t( x ' y; x '> y') exp(-D(x, y; xf, y')) 

where, we remark, amongst its several properties, the only relevant property of D for 
the present argument will be that 

D(x,y;x,y) = 0. (27) 

In the above thought-experiment, the result of any measurement of the statistical 
distribution of probe particles on the screen will, according to the naive pragmatic 
interpretation, arise from a formula of the form tr(/(y)p ut) where /(y) is the same 
function of the position of the probe particle that would be taken to represent the 
relevant observable in standard quantum mechanics. We have 

tr(/(y)Pout) = J J f(y)pout( x ,y, x ,y) dxd3 y 

But in view of fl26|) and (127]) this is the same thing as 

J J f(y)plut( x ^y^ x ^y) dxd3 y 

which is of course the same result which would be predicted for the same measurement 
in standard quantum mechanics. In particular, the statistical distribution of probe 
particles observed on the screen when we perform our experiment with the mirror in 
its ground state will be identical with that predicted according to standard quantum 
mechanics, which is what we set out to show. 
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One might of course choose to do a different sort of experiment to determine (p 2 ) of 
the state of our mirror - perhaps involving photons as probe particles and making use of 
optical devices such as interferometers etc. There is then, admittedly, some uncertainty 
as to what the predictions of the theory in [T] would be because of the relativistic 
nature of the photon and hence the inapplicability of the Newtonian approximation of 
[2]. But assuming (we shall return to this below) photons can, for this purpose, be 
treated as if they were non-relativistic particles (perhaps with a tiny mass) then e.g. an 
experimental arrangement which results in a detection of which of two optical channels 
a probe particle or photon goes down or an experimental arrangement which results in 
the examination of an interference pattern on a screen would, by these arguments, both 
give the identical result to that predicted by standard quantum mechanics, since both 
sorts of measurements are types of position measurement. 

We remark that similar considerations lead us to conclude that the results of 
experiments along the lines of those of Folman et al [8], to detect changes in parity of 
certain quantum mechanical states due to possible deviations from standard quantum 
mechanical laws, will, on the theory of [2] (and continuing to assume our reasoning 
continues to apply when some of the particles are photons) be identical to the results 
predicted by standard quantum mechanics, in spite of the fact that our naive pragmatic 
interpretation predicts changes in the expectation value of parity as in our formula (Tl8|) . 
The reason for this is again that the Folman et al experiment ultimately involves a type 
of position measurement (again of a photon). 

In fact quite generally, for any sort of experiment, as long as the result is obtained 
by a position measurement of some sort, to the extent that the approximation of [2] is 
valid, if we apply the theory of [2] to the closed total system consisting of the system of 
interest together with any relevant probe particles or 'pointer variables' etc. and adopt 
our naive pragmatic interpretation for this closed total system, then an adaptation of the 
above proof in the case of the bead example shows that the result of our measurement 
will be identical to the result predicted by standard quantum mechanics. To see this, 
suppose the total system consists of N bodies in all (protons, electrons etc.) and replace 
the variables (x, y) in the above proof by the iV positions of all of these. (It may be 
helpful to think of x as standing for, say, n 'system' variables Xi, . . . x n and y as standing 
for m = N — n 'pointer variables' yi, . . . , y m and we shall refer to this way of thinking in 
the next few parenthetical remarks.) Then the above proof goes through, thinking now 
of p? as the total would-be density operator at the time of preparation of the experiment 
and Pq U + as its Schrodinger time-evolute at the time of the measurement, /(y) as the 
appropriate position operator (i.e. an appropriate function of all the yi, ...y m ) an d 
taking the integrations to now be over the 3 N- dimensional total configuration space. 
The key point is again that the relevant decoherence exponent (now a function of the 
2N vector variables xi, . . . , x n ; yi, . . . , y m ; x' 1; . . . x^; y[, . . . , y' m ) will vanish on the 
diagonal (i.e. when x a = x' a Va 6 {l,...,n}, y& = y' b Vfe G {l,...,m}). This 
completes the proof of our general Position Measurement Theorem which we referred to 
in the Introduction. 
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In our corrected pragmatic interpretation we assume that a realistic measurement 
actually consists of a whole (Heisenberg-von Neumann) chain of measurements - each of 
which is a measurement in the sense we have been discussing up to now - as explained 
in the very well-known discussion in Chapter 6 of von Neumann's book [10]. Thus, for 
example, in the bead/mirror-probe-particle example with which we opened this section, 
the dynamics of the bead/mirror-probe-particle system are not the end of the story: the 
measurement of the position of the probe particle on the screen can in turn be analyzed 
by including the physics behind the formation of our spots on the screen. Let us suppose, 
for example that the spots consist of black grains of metallic silver formed when our 
probe particle hits a silver-iodide-coated screen. Then we would incorporate the screen 
too in the total quantum system and include suitable terms in the total Hamiltonian to 
describe the formation of the silver granules. Then one could perform yet another step 
of this sort to incorporate the physical mechanism by which light is used to measure the 
position of the silver-granule-spots and then by which the eye records the state of the 
light etc. etc. However far down such a chain one goes, our above Position Measurement 
Theorem will ensure that one will obtain the same result which is predicted by standard 
quantum mechanics, provided only that, at the stage of the chain at which one chooses 
to stop and apply the rule (Tl3|) . the measurement is a sort of position measurement. We 
could then simply declare < jxviiil > that the only realistic measurements are those where 
one does stop at a point of the von Neumann chain which consists of a sort of position 
measurement. Alternatively and again following a well-known line of argumentation (cf. 
e.g. [2H]) and cf. also the interesting recent paper [29J of Geoffrey Sewell) one could 
argue that, whatever may be the cause of the irreversibility, it is, in practice, only at a 
stage of the chain after something macroscopic and in practice irreversible has happened 
(such as the formation of our silver granules) that a measurement has really been made 
< jxix[ > and thereafter it is, in practice, immaterial whether the subsequent stages of 
the chain are analysed with quantum or with classical physics and, in view of that, it 
cannot make any difference if, insisting on still analysing quantum mechanically things 
which might be analysed classically, we decide to choose a stage of the chain in which 
the measurement is a sort of position measurement - e.g. a pointer reading on a dial. 

Our corrected pragmatic interpretation consists in adopting either of these points of 
view; it doesn't matter which. All that matters is that, for one reason or another, we take 
the view that a realistic measurement consists of a sort of position measurement at some 
point of the Heisenberg-von Neumann chain. With our above Position Measurement 
Theorem we then immediately conclude (as already stated in the introduction) 

On our corrected pragmatic interpretation, the experimental predictions of the 
Newtonian approximation of FJJ/ are identical with those of standard quantum mechanics! 

(A) 

Returning momentarily to our bead/mirror-probe-particle example, we remark 
that, actually, if the probe particle is sufficiently light (but we continue to use the non- 
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relativistic theory of [2]) then D(x, y; x'y') will be well-approximated by a function, say 
d, of x and x' only, so that 

PoutO, y; x ', y') - Pout( x > y; x '> y') ex p(-d(x; »'))■ 

Thus, for any observable A which refers to the probe particle so that, formally, it can 
be represented by a kernel A(y, y') which is only a function of y and y', we will have 

tr(^Pout)=y J J My^y')pont(x,y';x,y)dxd 3 yd 3 y' 
which by (j4]) is approximately equal to 

A (y> y')Pout(^) y'; x , y)dxd 3 yd 3 y' 
which, in view of d(x, x) = 0, is equal to 

^(y, y')Pout( x ' y'> x > y)dxd 3 yd 3 y' 
i.e. to 

= tr(Ap° ut ). 

This result will obviously generalize so as to conclude that, on our naive pragmatic 
interpretation, the results of any sort of measurement on a probe particle at some 
stage of a Heisenberg-von Neumann chain of measurement, will, if the probe particle 
is sufficiently light, be approximately the same as the results predicted by standard 
quantum mechanics. 

However, we wouldn't assign as much importance to this result as to our Position 
Measurement Theorem because it is only approximate and also because, in some given 
measurement chain, it may only apply e.g. at a stage of the measurement chain before 
classicality and irreversbility have set in. In such a case, we would still rather rely on 
our Position Measurement Theorem, applied at a later stage of the measurement chain 
where the measurement is, e.g. on a macroscopic pointer (which may well not be so 
light that one could make the above approximation) in order to justify the statement 
(A). 

As we have mentioned, all our results assume that the entire chain of measurements 
involves only non-relativistic particles, whereas, in practice (an example is the Penrose 
experiment we shall discuss in detail below) photons are often involved at some stage(s) 
of the chain. This is a gap which needs to be filled. However, we would expect that, 
provided (in the lab frame) the relevant photon states are not too energetic and provided 
they are sufficiently far from being on the verge of black-hole collapse then they may be 
treated as if they are non-relativistic particles with a tiny mass. In particular, we shall 
assume this to be the case for the photons in the Penrose experiment which we discuss 
next. 

Of particular interest to us are interferometer-like experiments along the lines 
of those suggested by Penrose [2TJ, [22] to detect the 'collapse of the wavepacket' 
of certain 'Schrodinger Cat-like' macroscopic quantum superpositions as predicted by 
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certain theoretical proposals for modifications of standard quantum mechanics, notably 
collapse-model proposals such as that of GRW [13] and the proposal of Penrose himself 
[20] concerning his gravity- induced 'quantum state reduction'. In one of its simpler 
versions <lxxl>. the setup proposed by Penrose is reproduced here in Figure 2. 




n 

Figure 2. Diagram of the Penrose experiment 

As described in |21j, a photon, emitted from a source, is passed through a beam- 
splitter and emerges in a superposition of two components with paths at an angle to 
one another, one, say 'component B\ of which is reflected directly back towards the 
beam splitter by a fixed mirror, whereas the other, say 'component A\ while eventually 
reflected back by a second fixed mirror, also reflects off a movable mirror of macroscopic 
dimensions and mass, both on its way towards and on its way back from that fixed 
mirror. The movable mirror is attached to an oscillator and (in a classical description) 
initially at rest. The mass of the movable mirror, the spring constant of the oscillator, 
the frequency of the photon and the spacing of the movable mirror and the second fixed 
mirror are tuned so that (staying with a classical description) 

in component 'A' of the superposition (in which the photon hits the movable mirror) 
its first encounter will set the movable mirror into motion, while its second encounter 
will occur at the movable mirror's initial location and bring it back to rest there - the 
movable mirror swinging to and from a macro scopically different location in between. 

(a) 

The positions of all the mirrors are also tuned so that, according to a standard 
quantum mechanical analysis, 

both components of the photon would arrive back at the beam splitter at the same time 
and with relative phases such that they reassemble into a state which travels back towards 
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the photon source, with no photons detected at the detector (see Figure 2) placed at the 
appropriate angle to that direction to catch any photons which might have emerged from 
the beam-splitter in its alternative direction. 

This standard quantum mechanical prediction depends however on the laws of 
quantum mechanics applying even though, at an intermediate stage of the experiment 
(i.e. between the two reflections of component A of the photon with the movable mirror) 
the total quantum state of the system is a superposition involving two macroscopically 
different positions of a macroscopic object, i.e. of the movable mirror. If a 'wave-function 
collapse' of this superposition occurs, e.g. as predicted by Penrose's 1996 proposal [20J, 
then one would predict that, in consequence, the coherence of the superposition of the 
two components of the photon beam as they reconverge on the beam-splitter will be 
lost and therefore that photons will arrive at the source and at the detector with equal 
probabilities. 

In this way, the Penrose experiment should be able to distinguish between standard 
quantum mechanics and the gravitational-induced quantum state reduction proposal of 
Penrose 1996 [20] (or between standard quantum mechanics and one or other collapse 
model). 

One might at first sight expect that Kay's 1998 theory (i.e. the Newtonian 
approximate theory [2] based on Kay's hypothesis [I J) would also predict a similar non- 
standard result for this Penrose experiment (for suitably large crystals etc.). After all, as 
we have seen, for suitably large crystals etc., that theory too predicts the decoherence of 
superpositions involving macroscopically different configurations. However, our general 
result above tells us (assuming, as we have discussed above, that the result still applies 
when the probe particles whose position is being measured are photons) that this can't 
be. Just as for any other experiment which ends in a position measurement, the 
prediction must be identical with the result predicted by standard quantum mechanics. 
This difference between the predictions of Penrose's 1996 proposal [20] and of Kay's 1998 
theory [2] for the Penrose experiment may be attributed to the following difference: In 
Penrose's work [20J it is assumed that once a macroscopic system enters a sufficiently 
macroscopic superposition (in the present case involving the movable mirror in two 
macroscopically different locations after the first reflection off it by Component A of 
the photon) then that superposition will collapse and remain collapsed irrespective of 
its future dynamics. (The same property holds in most collapse models.) In [2] on the 
other hand, the macroscopic superposition formed after the first reflection of component 
A of the photon at the movable mirror will decohere but then, at the second reflection 
recohere (but we emphasize that this is a feature of the Newtonian approximate theory 
[2], not of the full hypothesis of [1] see the remarks at the end of the following subsection 
and in Section [7]). We explain this point more fully in the following subsection. 
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4-1. Further discussion of the Penrose experiment 

In a mathematical description of the above standard quantum mechanical analysis, the 
time-evolving Schrodinger wave function of the photon/movable- mirror-cum-oscillator 
system may be regarded as an element of a total matter Hilbert space 

^matter = ^photon ® ^mirror- 

Suppose, in a wave packet description of the photon, the photon state at a time 
immediately after the photon emerges from the beam splitter takes the form -^=(7^(0) + 
7b(0)) G ^photon w here 7a(0), 7b (0) represent the A, B components respectively of the 
photon superposition at that time, and suppose that the initial state, /x(0) G 7^ m i rror , 
of the movable-mirror-cum-oscillator is its ground state, then we might interpret our 
statement (a) by the statement that 7^4 (0) <S> fi(0) G 7i ma tter wn l have evolved, by 
an intermediate time, t m (the suffix 'm' here stands for 'mid') shortly after the first 
reflection of the photon off the movable mirror, say, to a state 7a ® n(t m ) < dxxib > 
where fi(t m ), being a quantum description of the movable mirror in its macroscopically 
different location, will approximately satisfy (ji(t m )\fj,(0)) = 0, but that by the later time, 
tf (T here standing for 'final'), when the photon has undergone its second reflection off 
the movable mirror and is about to re-enter the beam splitter in the return direction, 
it will have further evolved to 74 (tf) <g> fi(tf) where, as the quantum counterpart to the 
classical statement that the movable mirror will have returned to its rest position, n(tf) 
will be taken to be equal to (a phase times) the ground state fi(0). On the other hand, 
Component B of the photon doesn't interact with the movable mirror at all, so the 
vector 7_b(0) ® /i(0) G T^total wn ^ evolve at any time into (a phase times) 7b(0 <8> A*(0) 
say. 

Putting all this together, standard quantum mechanics therefore would predict that 
the total state of the photon/movable- mirror-cum-oscillator system at time tf just before 
(both components of) the photon re-enter(s) the beam splitter will be < jxxiil > 

*Sltt d ef d (^/) = fy-r^f) + nttf)) ® MO) (28) 

so that the partial state of the photon itself (i.e. after tracing |^matt d r rd )^matter rd l 
over H m i rror ) will be the pure state with density operator 

Pphoton^*/) = l\Mf) + lB(tf)))((lA(t f ) + Mt f ))\ (29) 

and we assume that this is the photon state which passes entirely through the beam 
splitter in the direction of the source, leading to zero detection probability at the 
detector. 

Now, what we understand to be envisaged by Penrose [21], [22] < jxxiiil >. is that, on 
his 1996 proposal [20], the initial total density operator of the photon/movable- mirror- 
cum-oscillator system 

PWSaTter 6 = ^I(7a(0) + 7 b(0)))<( 7 *(0) + 7b(0))| ® \fl(0)) (/x(0) | 
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will (because of gravitational- induced quantum state reduction of the macroscopic 
movable-mirror-cum-oscillator superposition according to his proposal) evolve by time 
t m into a total density operator close to 

PmTttef^) = ^7A(O)(7A(OI®IMO)^(OI+^7fl(O)<7B(*m)l®l^(0))<//(0)|(30) 

and, because it is (at least implicitly) assumed that a reduced wave function cannot 
'dereduce' it is assumed that this will evolve by tf into 

Pmatter e (*/) = ^ Ita(*/)> <Ta(*/) | ® |p(0)><p(0)| + l -\ lB {tf))(l B {tf)\ ® |p(0))(/i(0)|. 

and therefore that the partial state of the photon itself as it re-enters the beam-splitter 
will be 

PpSon e (0 = ~l7A(t/)><7A(*/)| + \\lB{t f ))(l B {t f )\ 

which is expected to predict equal probabilities of detecting a photon at the source and 
at the detector. 

On the other hand, the theory of p], under the approximation of [2] would 
(see Notes <dvTb> and < lviiil >) require the replacement of the initial density operator 

l*St d e a r rd )(^St d ef d l> ^Stt^ as in ® on Matter = photon ® ^mirror by 
the partial trace, pj®^ of l^tal^total' over ^gravity where $^ is a suitable 
replacement for ^S t d e f d on H total = H matter ® ^gravity (= ^photon ® ^mirror ® 
gravity)- ^ n t ne formalism of [2] (cf. Notes <Jv3> and < Jviiit >) one would take (assuming 
as above that photons may be treated for this purpose as non-relativistic particles) 

*total^) to be given by 



= — ® M*) ® + 7b(*) ® MO) ® <?(0)) 



,Kay m _ J_ 
v/2 V 

where <?(t) is the (non-radiative) quantum counterpart to the classical Newtonian 
gravitational field of a (static) mirror at the position of the movable mirror at time 
t when the photon is in Component 'A' of its superposition. Platter (^) wm then clearly 
take the following values at times 0, t m and tf. At time zero, it will coincide with 
^matter rd (^)' at time tm ^ wm com cide with Platter 6 (^ m )' but at time */> ^ win a g ain 
coincide with fi^fe) = l*SSSte^(*/)><*SS*to d (*/)l ( see ©)• ^ other words, 
between times t m and tf, Patter wm ' rec °here' and in fact, as we know it must, will 
therefore result in identical experimental results to those predicted by standard quantum 
mechanics < lxxiv[ >. 

The reason for this recoherence is (cf. Note < jviiit > and also the last paragraph in 
Section E]) because the gravitational state-vector g(t) depends only on the state of the 
movable mirror at the time t and not on its history and therefore, since (reverting to 
the classical description which is relevant here) the movable mirror returns to its initial 
position at time tf, we will have g(tf) = g(0) < ixxvi >. Of course, in a fully relativistic 
treatment, this is expected no longer to be the case and in particular, if the motion of 
the movable mirror between times and tf causes the emission of one or more 'gravitons' 
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then one would expect (g(tf)\g(0)) = (approximately) in consequence of which the 
partial density operator p(tf) of the photon at time tf (and hence the experimental 
prediction) will coincide with Platter 6 ^/) - However, an order of magnitude estimate 
shows that, with typical experimental values for the mass and size and frequency of the 
movable mirror, the amplitude to emit a graviton will be tiny on time-scales relevant to 
the experiment < jxxvil > . 

5. A tentative interpretation in terms of events 

As discussed in the Introduction around the passage labelled (B) and discussed further in 
Notes < jxivl > . <|xyj>, < jxvil >. given any candidate for a fundamental alternative theory 
to standard quantum mechanics which is couched in terms of a time-evolving density 
operator, p(t), it would seem natural to attempt an interpretation in terms of 'events' 
which 'happen', by identifying the set of events which can happen at a given time t with 
the set of spectral subspaces of p(t) and identifying the probability with which a given 
such event happens with the associated eigenvalue, A, multiplied by its multiplicity, m. 
Thus, in particular, it is interesting to explore the prospects for such an interpretation 
in the case of the hypothesis of P] and one might hope that it would have a more 
fundamental status than either of our 'naive' or 'corrected' pragmatic interpretations 
(i.e. of the Newtonian approximate theory of [2]). From now on (and in Note < lxxviil >) 
we shall refer to it as our 'events' interpretation. To make a small beginning in this, we 
shall begin to explore this question here in the context of the Newtonian approximation 
of [2]. The caveats mentioned at the beginning of Section [L2l should of course continue 
to be borne in mind. (They will be discussed further Section [7]) . 

As we have already indicated in the Introduction and in the notes mentioned above, 
in the case of a would-be wave function for the centre of mass of a single ball in a 
Schrodinger Cat-like state 

where ipi and ip2 are sharply localized around two different locations, provided the mass 
of the ball is much larger than the Planck mass and the distance between the two 
locations a big enough fraction of the ball radius, the physical density operator will be 
well-approximated by ( fT2l) 

p= Idl 2 !^)^! + |c 2 | 2 |^ 2 )(^ 2 |. 

Our events interpretation then tells us that (except when |ci| 2 = 1/2 = |c 2 | 2 ) two events 
are possible, one corresponding to the subspace spanned by ipi, an d one corresponding 
to the subspace spanned by ?/> 2 . Further, these events will occur with probabilities |ci| 2 , 
|c 2 | 2 respectively. This is arguably just what one would hope for from a resolution to 
the Schrodinger Cat puzzle - see Note < lxiiil > for further discussion. 

In the remainder of this section, we shall mainly discuss how our events 
interpretation works out for a submodel of the one-dimensional Gaussian model which 
was introduced in Section CD, postponing discussion of more general models to Section 
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CD We recall from Section [T] that in this model the relationship between the would-be 
wave function, ip (on the line) and the physical density operator, p, is given by 

p(x, x') = ^{x)i)*{x')e- K{x - x ' )2 (31) 

and that the model is expected to give a good approximation to the physical density 
operator for the centre-of-mass motion of a uniform-mass- density bead of mass M and 
radius R constrained to move along a straight wire provided M is around or bigger than 
the Planck mass, or, in the case the wire terminates at hard stops a distance 2(R + 5) 
apart, for any mass provided 5 <^ R. 

We shall restrict ourselves further to would-be wave functions ip m this model which 
are either even or odd functions of x. As we show in the Appendix, using perturbation- 
theoretic methods, it is very easy to explicitly diagonalize any p in (I3T1) for a ip which 
satisfies this further restriction, to first order in k, and one finds, to this order, that 
there are only two non-zero eigenvalues (each non-degenerate), 

Ai = 1 - 2k(x 2 ), \ 2 = 2k(x 2 ) (32) 

with corresponding eigenfunctions 

01 = c l((l + K(x 2 ))lp — KX 2 ?p, (f)2 = C2X?p. (33) 

where C\ and c 2 are normalization constants. Here and below, (x 2 ) stands for the 
expectation value (ip\x 2 ip) of x 2 in the would-be state ip- 
As long as k(x 2 ) is very much less than 1, 

k(x 2 ) (= K(ip\x 2 ip)) < 1, (34) 

the eigenvalues (J32l) will be close to 1 and zero and therefore one expects this first-order 
perturbation theory result to give a good approximation to the true spectrum. 

We conclude that, if the would-be wave function, ip, satisfies (i/j\x 2 ip) <C 1 and 
is either an even function or an odd function of x then, to a good approximation, 
the possible events, according to our events interpretation, will be the one-dimensional 
subspaces spanned by (pi and 02 in (BSD anc ^ their respective probabilities will be Ai 
and A 2 in (132]) . To illustrate this, we have sketched, in Figure 3, 4>i and (f> 2 when the 
would-be wave function is the ground state ( |24l) 

tjj( x ) = 5- 1/2 cos(irx/25) 

of our bead-on-a-wire example (see Section [1] after equation <Q and equation (124|) in 
Section H] and below) with stops at ±(R + 5). 

In fact, as explained in Section [T], as long as 5 R then the Gaussian model will give 
a good approximation for the physical density operator, p, of this bead-on-a-wire-with- 
stops example whatever the mass, while, in view of (I35D and, recalling that k = 9M 2 /R 2 , 
we expect the above first-order perturbation theory results to be applicable provided 
5 R/3M which is of course not an additional restriction as long as M is around or 
smaller than the Planck mass. 
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Figure 3. A typical wave function for the one-dimensional Gaussian model and its 
two 'events' 



On the other hand, we note that for the application to our bead moving along an 
infinite straight wire, then the condition fl34l) amounts, in view of k = 9M 2 /R 2 , to the 
condition 

(iplx 2 ^ < R/3M. (35) 

so, since, as we saw in Section [[J we need to assume M to be around or greater than 1 
for the Gaussian model to be a good approximation at all, we can trust our above first- 
order perturbation theory result only when the width of the would-be wave function for 
the centre-of-mass of the bead is very much less than the bead's radius. 

In any case (i.e. whether or not the one- dimensional Gaussian model is a good 
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approximation for some bead-on-a-wire system) it is interesting to notice (see Figure 
3) that (to first order in k) if we take, say, an even would-be wave function if>, then 
0i is even and approximately resembles if) and occurs with a probability, Ai, close to 
1, but there is a small but non-zero probability, A 2 , for 2 to occur and this is a quite 
different function from (in fact it is, of course, orthogonal to) ip and it is odd! Moreover, 
the expectation value of parity (V), by which we now mean the ordinary appropriately 
weighted sum of probabilities (+l)Ai + (— 1)A 2 is then, by fl32|) . given by 

(p) = 1 - 4k(^|x 2 ^>. (36) 

and coincides with the formula (1201) for the expectation value of parity (albeit defined 
there, differently, to be tr(pV)\) obtained under the naive pragmatic interpretation in 
Section [21 (And one can similarly reproduce (12TT) when ip is odd.) The coincidence is of 
course explained by the fact that, for a would-be wave function, if), which is either even 
or odd, the conventional Parity operator V is simultaneously diagonalizable with p and 
has eigenvalue +1 on 0i and —1 on 2 < lxxviil >. One can also show, using the results in 
Subsection I9.2l in the Appendix, that similar coincidences obtain, for the analogous pair 
of different possible meanings for the phrase 'expectation value of parity' in the case of 
the three-dimensional situation discussed in the sentence containining equation [18] and 
the subsequent sentence. 

6. Entropy and its time dependence 

As soon as one knows how to diagonalize p, one can calculate its (von Neumann) entropy, 
S = — tr(plnp), and other entropy-like measures of its amount of decoherence - in 
particular Si = 1 — tr p 2 (see e.g. [21]) which, while less significant physically, is often 
more tractable mathematically. 

Staying with the one-dimensional Gaussian model and working to lowest order in 
k as explained in Section [5] and the Appendix, one finds, for a given would-be wave 
function if), 

S = -2K(x 2 )\n(2eK(x 2 }), 
Si = 4k(x 2 ) 

where, as in the previous section, (x 2 ) stands for (if)\x 2 if)). These are expected to be 
good approximations as long as k(x 2 ) is very much less than 1. 

If the would-be wave function if) evolves in time according to some would-be 
Hamiltonian H, as discussed in the Introduction around equation ([8]), so that the 
physical density operator evolves according to the master equation ([TO]) 

p = -i[H,p] - -£[x, \p,p\], 

then, clearly, the time-evolution of the approximate expressions for 5* and Si above 
is controlled by the time-function t i— > (if)(t)\x 2 if)(t)) . (Below, we call this (x 2 (t)).) 
If this increases, our approximate S% (which is just 4k times this quantity) will 
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obviously increase and so will our approximate S provided k{x 2 ) < 1/e 2 . (But we are 
anyway assuming that k(x 2 ) < 1.) If ip(t) evolves according to the free Hamiltonian 
H = p 2 /2M, we have the easily derived and well-known 'spreading-of-the-wavepacket' 
result that (x 2 (t)) will take the quadratic form 

(x 2 (t)) = A + (2E/M)(t-t min ) 2 (37) 

where E is the (conserved) expectation value (ip(t)\Hif)(t)) of the energy in the state 
ip{t). The feature of ( 1371) to which we wish to draw attention here is that (x 2 (t)) has a 
single minimum (with value A) - occurring consequence of this, substituting 

( 1371) into either of our above approximate formulae for S or Si gives rise to a 'two-sided' 
entropy- increase result in the sense that, if the time t m j n is taken as an appropriate 'zero 
time' then the entropy will only increase on either side of that (in other words, it will 
decrease monotonically before the zero time and increase monotonically after it) at least 
until K,(x 2 (t)) ceases to be very much less than 1, at which point our approximations 
and assumptions will break down. 

This two-sided entropy increase result would be applicable, for example, to the 
centre-of-mass degree of freedom of our bead of mass M (see Section [5]) bigger than or 
around the Planck mass, free to slide along an infinite wire and with an initial wave 
function at t m i n such that A in ( 1371) is very much less than R 2 /9M 2 , where R is the 
bead radius, as long as the times considered are sufficiently short for the wave-packet 
not to have spread so much that the condition (135]) ceases to hold. To give a quantitative 
example, consider a Planck mass bead 2 x 1CT 5 g) with radius R = 10~ 2 cm and a 
state in which (a; 2 (t min )) = A = KT 9 cm 2 . Then S lmin (= An A = Ax^-A) « 4x 10" 4 . 
Moreover, if we assume that we can trust ( 1321) (and hence our formula for Si) as long 
as (x 2 ) < 10~ 7 cm 2 , say, so that k(x 2 ) < 1CT 2 , then the time, t, at which (x 2 (t)) attains 
this value, can be estimated (by setting ^r^f^ 2 = 10~ 2 and taking E = ^ with p 2 
by an uncertainty principle estimate = jt) to be, within an order of magnitude or so, 
around 10 9 years! ... by which time, Si will have increased to around 4 x 1CT 2 ! 

It is interesting to compare and contrast this sort of (two-sided) entropy increase 
result (for either S or Si) with the sort of entropy increase results which one can obtain 
for the BLP master equation [27] 

P = -i[H, p] - c[x, [x, p}} (38) 

(Here H is some choice of Hamiltonian and c a positive constant.) which may be regarded 
as a prototype example of a master equation of GKS/L form. (See Note <Hxb>.) 

For this BLP model, it is easy to show that, irrespective of the Hamiltonian H, the 
quantity Si will, for positive times, increase monotonically in time. To show this, begin 
by noticing that it follows immediately from the definition Si — 1 — tr(p 2 ) that the time 
derivative, Si = — 2trpp = — 2tr(— i[H, p]p) + 2ctr([x, [x,p]])p). 

The first term here vanishes by the cyclic invariance of the trace. The second term 
equals — 2ctr([x, p] 2 ) = +2ctr([xp}^[xp}) which is > 0. (End of proof.) 

This situation differs in at least two important respects from that we saw above 
for our Gaussian model, which, we recall, satisfies the master equation (ITOl) when H is 



Quantum gravitationally decohered quantum mechanics 



30 



taken to be the free Hamiltonian H = p 2 /2M. To explain the first difference, we first 
need to notice that the BLP master equation fl38|) is itself of a quite different nature to 
( fTOl) in that (see Note <drxb>) it holds unrestrictedly for any p but only for positive times, 
whereas (TTOT) is only meaningful for p for which p(x, x') takes the special form ([3]) but 
then holds both for positive and negative times. (See Note <E0>.) This difference is of 
course reflected in the different sorts of entropy increase results which are appropriate, 
with our two-sided result for ffTOl and the above, more familiar, one-sided result for fl38|) . 

Second, the above 'one-sided' entropy-increase result for the BLP master equation 
( I38|) holds irrespective of the choice of Hamiltonian H in (138|) . In contrast, the two- 
sided entropy-increase result for our master equation (ffOl) which we obtained above 
in the case that H in f fTUj) is the free Hamiltonian H = p 2 /2M, depended on the 
form of the Hamiltonian. Indeed, suppose one were to take, instead, for example, the 
Harmonic oscillator Hamiltonian H = p 2 /2M + kx 2 /2 (describing, say, our bead-on-a- 
wire when the bead is connected to the origin with a spring with spring constant k) 
then it is obvious that, since the would-be wave function ip(t) will be periodic in time 
(with period 27r(M//c) 1//2 ) then so will Si and S - these being functionals of ip. This 
second difference can be traced (see Notes <Hxb> and <H>) to the difference that master 
equations of GKS/L form such as BLP arise from tracing over 'radiation' type modes 
which 'fly away' whereas our master equation (flQj) arises (cf. the last paragraph of 
Subsection 14.11 and Notes < jviiil > and < ixxvi >) from tracing over modes which (in our 
bead-on-wire interpretation) are dragged around by the would-be wave function of our 
bead. 

The relevance of the two-sided entropy result obtained in this Section for our 
Gaussian model to the general understanding of the Second Law entailed by the 
hypothesis of pQ will be discussed in the next (Discussion) section. 

7. Discussion 

Our main purpose has been to explore further the non- relativist ic approximate theory of 
[2] which, as argued in [2] and explained in the Introduction, is expected to be the non- 
relativistic limit of any full theory of quantum gravity consistent with the hypothesis of 
p] which, in its turn, as argued in [1] (see also Note <EJ>) offers a natural resolution 
to the several puzzles and problems mentioned at the start of the Introduction and in 
Note <S>. 

As we have explained in the Introduction, this approximate theory amounts, in 
the most general form in which we have stated it (see Equation ([7]), Note <jviit> . and 
the parenthetical remark after Equation ([8])) to the following modification of standard 
(many-body) non-relativistic quantum mechanics: 

One has an 'underlying' many-body wave function \jj(t;xi, . . . , xjv) which evolves 
in time exactly as in standard non-relativistic quantum mechanics for the usual many- 
body Hamiltonian that one would adopt for the problem of interest. However, while 
one would normally take the density operator of the system at time t to be given by the 
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(pure) projector onto \l/, which in its position-space representation is given by 

p (t; xi, . . . , Xat; x'i, . . . , xV) = Xi, • • • , xjv))^(t; x'i, . . . , xV)* (39) 

one declares that the physically relevant density operator p is given, in its position-space 
representation by 

pit; x 1; . . . , X7V ; xi, . . . , x'J = p (t; x ls . . . , x w ; x^, . . . , x ^) e -J>(*ir».**;*i,...,*W (40) 

where the decoherence exponent -D(xi, . . . , x^; x' 1; . . . , x^) is a specific function on 
configuration space (see Note <lviil>) which depends both on the mass of each body and 
also on a suitable radius (see the Introduction, Note <dvnl> and [3]) which is ascribed to 
each particle and which satisfies the general properties stated after Equation <^7§ - the 
most important of which for the present discussion are 

(a) D tends to be large for pairs of configurations (xi, . . . , xjv), (x'i, . . . , xV) which 
(in the terminology of [2]) involve 'mass relocations' around or bigger than the Planck 
mass. As discussed further in [2], for a given pair of configurations, the mass relocation 
involved may be defined to be the minimum amount of mass that has to be moved 
in order to convert the unprimed configuration into the primed configuration (or vice 
versa) . 

(b) D vanishes on the diagonal, i.e. 

.D(xi, . . . ,xjv;xi, . . . ,Xjv) = 

As was discussed in the Introduction and in [2], Property (a) leads to the 
suppression of macroscopic ('Schrodinger Cat-like') superpositions while property (b) is 
what lies behind our Position Measurement Theorem which, on our corrected pragmatic 
interpretation (see Section Hj) leads us to conclude that the experimental predictions 
of our non-relativistic approximate theory are identical to those of standard quantum 
mechanics. (As we cautioned at the beginning of Subsection 11.21 this conclusion is also 
predicated on being able to treat the relevant small non-relativistic systems as closed 
for the purposes of this discussion of measurement. We shall provisionally assume that 
it is and give arguments which confirm this towards the end of this section.) 

As we remarked in the Introduction (see before the paragraph labelled (A) in Section 
II. 2p the latter conclusion might seem surprising. After all, for large macroscopic systems, 
the difference between the mathematical quantities p and po can be huge by any measure. 
On the other hand, the conclusions of [2] that quantum gravitational effects could have 
any sort of big effect on every-day laboratory sized quantum systems were themselves 
surprising. So it is perhaps reassuring that we have now arrived at the conclusion that 
(in the non-relativistic regime) they don't have any measurable effect after all. The way 
in which they manage not to have a measurable effect, however, is rather subtle and 
despite the experimental indistinguishability, the theory of [2] possesses properties and 
prospects quite different from those of standard quantum mechanics as we now discuss. 

First and foremost, the theory of [2] achieves the goal of replacing standard 
quantum mechanics by a modified theory in which macroscopic ('Schrodinger Cat-like') 
superpositions are absent. As far as the achievement of this goal is concerned, it would, 
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indeed, seem to be an attractive rival e.g. to the 'collapse models' such as that of 
GRW [13] and others, which achieve the same goal at the expense of certain ad hoc 
modifications of quantum mechanics (and which also entail certain small modifications 
to the experimental consequences). 

It is interesting to notice that, in the light of the present work, one realizes that 
it would have been possible to achieve the same goal simply by positing that the usual 
pure density operator fl39|) of standard quantum mechanics should be replaced by a 
physical density operator given by ( )40|) with D any (ad hoc) choice of function on 
Configuraton space x Configuration space (perhaps quite unrelated to quantum gravity) 
designed to have the properties (a) and (b) above (and the other general properties 
stated after Equation ([7])). For, once one posits any such replacement, one will be 
able to prove a Position Measurement Theorem and to propose a corrected pragmatic 
interpretation etc. To the best of our knowledge, however, no such proposal had hitherto 
been made. In any case, the proposal with our approximate theory of [2] would (in the 
absence of any current experimental guidance) seem to be more satisfactory than that 
with any such ad hoc choice of D and also more satisfactory than the resolution achieved 
with any ad hoc collapse model, not only because it is less ad hoc, but also insofar as 
it is part and parcel of the hypothesis of pQ which offers a much more wide-ranging 
collection of simultaneous resolutions to several other puzzles as explained at the outset 
of the Introduction and in Notes <E>, <dmt>. <Hvt>. <EJ> and <jxiib>. 

The fact that our non-relativistic approximate theory manages to resolve the 
Schrodinger Cat puzzle (i.e. by eliminating all macroscopic Schrodinger Cat-like 
superpositions) while preserving (on our corrected pragmatic interpretation) exactly 
the same experimental predictions as standard quantum mechanics might be regarded 
as another advantage. But it also of course means that, viewed in its own right, it is 
only what one might call a philosophical resolution in that it leads to a changed ontology 
(i.e. with no Schrodinger Cat-like macroscopic superpositions) from that of standard 
quantum mechanics but makes no experimentally distinguishable predictions. 

However, it is saved from being 'merely' philosophical since the experimental 
indistinguishability of course only holds for our non-relativistic approximate theory. 
As is implicit in Section H] (see also Note < lxxvil >) when post-Newtonian corrections 
to our theory become significant, one would expect experiments along the lines of the 
Penrose experiment to detect differences. We next turn to consider at more length what 
the post-Newtonian corrections to our theory might be like. 

As our starting point, we return to the Newtonian time-evolution rule summarized 
in the first paragraphs of this section. When some details are stripped away, this 
specifies, for a given zero-time density operator p(0) with position space representation 
as in (|39|) with t set to zero, a map 

1 1 ► pit) (41) 

from the full real time-line (i.e. the full set of positive and zero and negative times) to 
the space of density operators on the relevant many body Hilbert space ^matter- 
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This rule will arise as the solution to a suitable master equation. We won't write 
this down explicitly; it is straightforward to obtain it by differentiating Equation (1391) 
with respect to time just as ([9]) is obtained by differentiating the one-body special case 
(JHJ). Instead, we will write it schematically as 

p = -i[H, p] + Newton p (42) 

where Newton is the appropriate operator which acts on the set of density operators. It 
will be a 'master equation of MDM type' in the terminology introduced in Note <KJ>. 

In the post-Newtonian regime, we would expect the nature of the dynamics of the 
physical density operator p(t) (i.e. of the trace over 7^g rav ity °f the time-evolving pure 
total density operator as explained in Subsection II. ip to be qualitatively different (as we 
are about to discuss) but we would still expect it to admit of a description in terms of 
a map of the form of (|4ip - i.e. from the full real time-line to the same space of many- 
body density operators. However, unlike ( )40l) and its associated master equation f j42l) . 
one expects a new sort of rule determining the time-evolution of p(t) which incorporates 
the back-reaction on our many-body matter from radiated gravitons. We would not 
necessarily expect this new rule to obey an exact master equation but, in suitable 
models, we would expect there to be a preferred zero-moment of time at which the 
entropy is a minimum - the many-body post-Newtonian counterpart to the moment 
(^min m (0Z])) of minimum entropy in our one-dimensional Gaussian model example 
of Section [6] - and, while, of course, one now expects the rule which determines p(t) 
at other times to depend on the underlying total time-zero state and not just on the 
physically relevant density operator p(0) (which is its partial trace over gravity - see the 
beginning of Subsection 11.11) . we would expect, in view of the considerations discussed 
in Note <Jixt>. that, for positive times p(t) will approximately obey a master equation 
of schematic form 

p = —i[H, p] + Newton p + postNewton p (43) 

where postNewton is approximately of GKS/L form while, for negative times it will 
approximately obey a similar master equation but with the term postNewton replaced 
by a term antipostNewton which is of anti-GKS/L form <dixb>. The reason is that (for 
positive times) the post Newtonian terms will arise by, at each time, tracing the full 
matter-gravity (pure) density operator over modes of the gravitational field which are 
radiative (see <dixl>). and, unlike the non-radiative modes which, in virtue of being 
'dragged around' by the matter are responsible for the special features of Newton as we 
discussed in Subsection 14.11 Section [6] and Note <G0>, radiative modes will, of course, 
be expected to fly away from the matter. (And there is of course a similar argument 
for the expected anti-GKS/L property for negative times.) Moreover, while the terms 
postNewton and antipostNewton are expected to be negligibly tiny in comparison to 
Newton in the Newtonian regime (see Note < jxxvil >). we would expect that, as one 
moves away from the Newtonian regime, they will become comparable in size, and as 
one approaches an extreme relativistic regime, they will be much more important. So 
one might say that the non-relativistic approximate theory, which is the main subject 
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of the present paper, is expected to just be the small non-relativistic tip of a large 
relativistic iceberg! (Of course, the possibility of describing p as a many-particle density 
operator will presumably break down at some point and need to be replaced e.g. by a 
description in terms of quantum fields or strings etc.) 

We next turn to discuss the issue of which systems deserve to be considered closed 
for our various purposes. As we have indicated towards the beginning of Subsection 11.21 
and warned in several places above, it was already concluded in [2] that, at least as far 
as the concept of entropy is concerned, ordinary small laboratory-sized non-relativistic 
systems cannot be considered closed. We shall next recall the reasoning behind this 
conclusion and shall then discuss whether or not such a conclusion should also apply 
to the other matters under discussion here - notably our conclusions about the Penrose 
(and other) experiment (s) and our tentative interepretation in terms of events. 

By a well-known discussion (see e.g. the article by Joos in |26j) based on the 
argument that there will always inevitably be quantum mechanical entanglement over 
widely spatially separated subsystems, one can argue that it only really can make exact 
sense to regard a subsystem of the universe as closed and to regard its underlying state 
as describable by a pure (i.e. vector) state if the subsystem is the entire universe. 
Indeed, at a classical level, it was famously pointed out by Emil Borel in 1914 that 
the motion of a few grams of material on the star, Sirius, would, just because of its 
gravitational interaction, after a short time, significantly change the configuration of 
a gas in a box on the Earth. Similar considerations at a quantum level (and for not- 
necessarily gravitational interactions) make it plausible that, in a quantum description, 
in the actual state of the universe there will be a considerable degree of entanglement 
between a bit of gas in a box on the Earth and a bit of gas in the star Sirius etc. (For 
a discussion of all this, see for example Chapter 3 by E. Joos in |26j) 

However, one can ask whether it is legitimate to assume a given system is closed for 
the purposes of some given calcuation - if one is content with a good approximate (rather 
than an exact) answer. We know from experience that, for many questions, such an 
assumption can be successful for laboratory-sized systems. Otherwise standard quantum 
mechanical calculations would not be anything like as successful as they actually are. 
However, if our purpose is the explanation, on the hypotheses of [lj, and within the 
approximate Newtonian framework of [2j of how e.g. a box of gas in the laboratory 
can have an entropy comparable to the usual entropy one would calculate for a thermal 
equilibrium state at say room temperature, then it appears, from the preliminary model 
calculations reported in j2] that this cannot be done: The typical entropies we calculate 
(see [2]), while non-zero, are orders of magnitude smaller. Our conclusion from this 
is not that the theory of [2] or the calculations are wrong, but rather that hot boxes 
of gas in the laboratory cannot be regarded as closed systems. Rather, as we discuss 
further in Note <jxiib> (see also Figures 4 and 5) we conclude that the entropy of such 
systems needs to be explained by regarding them to be open systems and is due to 
tracing over a total environment, part of which is gravitational, but the main part of 
which is essentially the same matter environment which (on one standard view about 
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the origin of the entropy of gases in boxes - see the last paragraph of Note <lxii[>) would 
usually be held responsible for their entropy. 

It then becomes an important question, whether, to get physically correct values 
for the entropy of such laboratory-sized systems and to get physical thermodynamic 
behaviour (i.e. a 'Second-Law' result to the effect that the entropy of a closed system 
always increases) it could suffice to regard as closed, say, the solar system, or whether 
one needs to include our whole galaxy or whether one actually needs to go as far as the 
entire universe. We have argued in [2] and will argue further below that any answer 
to this question must involve enlargement to a system, parts of which are essentially 
relativistic in nature (the remark in Note < jxiib > about total closed sytems which contain 
black holes is also relevant to this conclusion). Thus the galaxy or the entire universe 
would seem to be the more likely answers. And if the galaxy turns out to be a valid 
answer, then one might speculate that if it is true, as is currently believed likely, that 
every large galaxy contains a supermassive black hole, then this fact may play a role in 
this. 

What we have studied in [2] and in the present paper are small laboratory- 
sized model systems on the (we have concluded above as far as entropy is concerned) 
counterfactual assumption that they are closed. While not expecting their entropies to 
be physically realistic, we can endeavour to learn lessons about the mechanisms which 
give rise to entropy and about the way that entropy tends to vary with time. In [2] 
we considered a number of laboratory-sized model static systems and calculated their 
entropies on this assumption, and showed, in particular, that these entropies, while 
(as we mentioned above) small, appear to correlate, to some extent, with their degree 
of 'matter clumping'. In Section [6] here, we have augmented these results with some 
first results on the time- dependence of entropy and found, for our Gaussian model to 
first order in k (describing our free bead-on-a-wire model) a two-sided entropy increase 
result whereby there is a moment of minimum entropy while, as a result of 'spreading 
of the would-be wave-packet', entropy increases for times on either side of this moment 
(i.e. it decreases before this moment and increases after it). For more 'realistic' small 
non-relativistic many-body systems we would expect there to be a similar two-sided 
entropy- increase result in suitable circumstances, but we expect the mechanism behind 
it to involve not only spreading of the would-be wave packet, but also, and in fact 
often mainly, 'dynamical clumping', by which we mean changes in time of the amount 
of matter clumping. Turning to the post-Newtonian regime, we would still expect to 
obtain a two-sided entropy-increase result in suitable models, but we expect the way 
one will understand mathematically the origin of this result will be different (and this 
will be traceable to the different behaviour of the non-radiative and radiative modes 
which will be responsible for the two regimes - see <Hxb> and <|x]>). In the case 
of our Gaussian model and (one assumes) in the case of non-relativistic many-body 
systems, the demonstration of such a two-sided entropy increase result depended/will 
presumably depend on a study of the relevant single exact master equation of 'MDM' 
type (see after (142!) which describes the dynamics of p{t). In the post-Newtonian regime, 
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one expects, instead, that after choosing a suitable model for a low-entropy inital p(0), 
the demonstration of a two-sided entropy increase result will depend on a study of 
an appropriate approximate master equation of form ( )43l) for positive times and its 
negative-time counterpart for negative times. To the extent that one can neglect the 
Newton terms in these (cf. our 'tip of iceberg' remark in the paragraph containing ( )43l) ) 
the question will be whether or not the GKS/L term postNewton satisfies the criterion 
of Benatti-Narnhofer (36] discussed in Note <drxb> (etc. for negative times). 

The interpretation of such two-sided entropy increase results will of course be that, 
while, theoretically, they represent possible dynamics for all times, positive and negative, 
in actual physical application, only the positive-time behaviour will be relevant and the 
negative time part should simply be discarded. Either one has in mind that one's closed 
system is modelling the whole universe in which case p(0) will model the initial state of 
the universe, or, if one is interested in a model which idealizes some subsystem of the 
universe as closed (e.g. a star collapsing to a black hole in an otherwise empty universe 
as is done in Note <EJ>) then p(0) will represent the initial state of that system and 
its low entropy can either be thought of as an assumption of the model or alternatively, 
perhaps, as ultimately explainable outside the model as traceable back, via processes in 
the wider universe, to the low entropy of the initial state of the universe. (Cf. the notion 
of 'branch systems' in [30] although note that, in view of our remarks above about how 
large closed systems may have to be, it is a delicate question if/when/how this notion 
is applicable here.) 

Of course, the only result we have proven here is our two-sided entropy result, 
obtained in Section [6j for our one- dimensional Gaussian model, interpretable in terms 
of a model closed system consisting of our free non- relativist ic bead on a wire. This is of 
course a thoroughly counterfactual model closed system but nevertheless our result is of 
value as a first entropy-increase result on the hypothesis of pQ . If that hypothesis and the 
extra assumptions involved in [2] are correct, then that model will, as we have indicated 
above, be of some significance as one extreme tip (it concerns only one single small 
non-relativistic bead) of a non-relativistic many-body tip (cf. the paragraph containing 
( 143]) ) of a full relativistic iceberg which is the Second Law of Thermodynamics! 

The other side of the coin to entropy increase (on the hypothesis of p]) is of course 
dynamical decoherence since, after all, aside from its role in thermodynamics, our S can 
be thought of as a measure of the 'amount of decoherence' of our closed system. So, 
if our expectations concerning results of Second-Law type for suitable (say, laboratory- 
sized non-relativistic) models are fulfilled, then we will be able to say that, for suitable 
'time-zero' initial total states, the amount of decoherence will always increase. Thus, 
returning to our comparison with theories of 'collapse model' type, we will be able 
to say that our non-relativistic approximate theory shares with 'collapse models' (i.e. 
in addition to eliminating Schrodinger Cat-like superpositions as discussed above) the 
ability to predict a continual process of decoherence. (The way that GRW does this 
is mentioned in Note <drxb>.) However, the comparison of the two sorts of theories has 
to be more complicated than that, since, as we have argued above, small laboratory 
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systems cannot, on our theory be regarded, for the purposes of calculating entropy, as 
closed. Regarding them, instead, as open (as discussed in Note <dxiH>) we see that, 
instead, our theory (i.e. of [2]) entails, as far as its predictions of decoherence are 
concerned, essentially the same predictions as (see Note <EJ>) the usual 'environment 
induced decoherence paradigm'. (See however the next-but-two paragraph.) 

Turning to our tentative notion of 'events' defined according to the passage labelled 
(B) in Section (II. 2p . we have to expect that, for related reasons to those given above, 
the detailed structure of the set of events and their probabilities and of the way these 
change with time will differ considerably, for a large relativistic closed system, from 
the structure etc. of this set etc. when one analyses simple non-relativistic models 
of closed systems as in Section O Nevertheless, the results of Section [5] are still 
of interest as a model within which one can attempt to explore what might be the 
mechanism by which such an events-based interpretation might be able to supplant the 
stop-gap naive and corrected pragmatic interpretations discussed in Sections [2] and HI 
What one hopes for here is to have a theory which explains how quantities such as 
position and momentum and parity etc. which are treated as 'observables' in our naive 
and corrected pragmatic interpretations emerge in some suitable sense from an events 
interpretation. Our discussion, in Section [5] (see also Note < lxxviil >) of the relationship 
between equations (1201 and (136|) in the context of our one- dimensional Gaussian model, 
appears to give a promising indication, for the example of 'parity' of how this can 
happen and the challenge which remains is to study whether an explanation along the 
lines indicated by that example, of the emergence of less primitive observables in more 
complicated (many-body) examples. 

We emphasize that even if a theory along such lines for how observables emerge 
from our events interpretation is possible, it is not obvious that it will be possible 
staying within the context of small non-relativistic model closed systems. After all, 
we argued above that for questions related to entropy and thermodynamics, such 
models need to be rejected in favour of models of large relativistic closed systems. 
However the small success of our discussion of parity encourages us to hope that, 
for these different purposes, our small non-relativistic models may suffice. Indeed, it 
seems conceivable that, even though a different, and perhaps richer and ultimately 
arguably more 'physical', theory of events would be had for large relativistic models 
of closed systems and by treating actual ordinary laboratory-sized systems as open 
subsystems of these, the theory obtained by modelling actual ordinary (non-relativistic) 
laboratory-sized systems as closed may, in itself, yield a viable events-based replacement 
for standard quantum mechanics. 

Concomitantly, one might possibly argue that, returning to our comparison between 
our theory and 'collapse models', if we treat a given small non-relativistic laboratory- 
sized system as if it is 'closed' for the purposes of calculating its entropy, 5*, then, 
while the resulting S should not be taken seriously as having anything to do with 
thermodynamics or 'true extent of decoherence', it might still be an interesting measure 
of 'decoherence relevant to the model' and might, under certain circumstances, always 
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increase. 



Finally, returning to our naive and pragmatic interpretations which we discussed at 
the outset, we return to the question whether it is justified, for the results which were 
obtained within them (including our Position Measurement Theorem and including our 
analysis of the Penrose experiment) to be taken to be physically relevant. The question 
is whether such matters can (unlike what we argued above for questions of entropy and 
thermodynamics) be discussed in the context of the models of small non-relativistic 
closed systems which we have considered here. Another way of putting this is to ask 
whether, e.g. in our discussion of the Penrose experiment, we can ignore entanglement 
of our small model closed systems with the environment. The answer would seem to be 
that it is justified since (cf. Note < lxiil >) the issue of entanglement of such systems with 
the environment within our theory is more or less the same as what would usually be 
regarded as the issue of the entanglement of such systems with the environment and this 
is, of course (see especially [9]) precisely the issue which makes the Penrose experiment 
so difficult (but not impossible) to perform and which the experiment is designed to 
overcome. 

(A brief summary-statement of some of the main conclusions from the above 
discussion with an emphasis on comparing and contrasting our Newtonian approximate 
theory with 'collapse models' such as GRW [13] was given towards the end of the 
Introduction.) 
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9. Appendix: Approximate diagonalisation of the density operator in the 
one-dimensional Gaussian and some other, related, Models 

Given p as in ( {TBI 



p{x,x') 



where 



po(x,x') = 4>(x)ip*(x'), 
our aim is to find the eigenvectors and eigenvalues of p. I.e. we wish to solve 

P |0) = A|0), 



(44) 



Quantum gravitationally decohered quantum mechanics 39 

for |0) and A. This could be studied exactly, but we content ourselves here with solving 
the simpler problem where p is replaced by its first order in k form, which shall call p K , 

p K (x, x') = po(x, x') (l — k(x — x') 2 ) = po(x, x') (l — kx 2 + 2kxx' — kx' 2 ^ 

= p (x,x') + Kp!(x,x'), 

where 

Pl(x,x') = —Klp(x)lp*(x')(x — x') 2 

and the resulting eigenvalue equation is solved to first order in k. It is useful to note, 
that, in Dirac notation, p K may be written 

p ~ p + K px = \ip)(ip\ + k (-x 2 \^)(ip\ + 2x\il)){ip\x - \ip)(ip\x 2 ) . (45) 

Taking this as our starting point, we can make the following perturbation-theoretic 
argument: If k is small, we expect p, and hence also p K to have an eigenvector close to 
and that its eigenvalue will be close to unity Thus we write this first eigenvector 

as 

|0 1 ) = |^) + /c|V±>, (46) 
where \<f>i) is not yet normalized and we assume the orthogonality condition, 

holds. As tr p = tr p K = 1, and p is positive (we will confirm that p K is positive for 
sufficiently small k below) we expect that the eigenvalue of \<pi) will be slightly less than 
unity so we write 

p"\<f> 1 ) = (l-<wc)|^>, (47) 
where a is to be determined. From (145|) . (|46l) . and (1471) we get 

P K i0!) = (p + m) m + = (i - o«) m + k\^)) . (48) 

We can obtain |-0j_) from (T481) by expanding out the brackets to give 
|V>±) = (a + pi)|V>), 

and substituting the expression for pi, given in ( )45l) . into the above equation gives 

|^ ± ) = a \ip) - x 2 \ip) + 2(x)x\if>) - (x 2 )\ip). (49) 
By only keeping terms up to order k in (1481) we get 

rcpil^) = -arel^) + ^IV'-l) (50) 
Now acting on (1501 with it can be seen that 

a = -(^|pi|^)=2(Ax) 2 , (51) 

where 

Ax = ((V>k 2 |^> - <-0l^l^> 2 ) ^ 



■2 
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Note that the above is the uncertainty in position of the would-be wave function (as we 
call it in Subsection II. ip . So we are now in a position to write an explicit expression 
(not normalized) for our first eigenvector of p to first order in k, using f T46l) . ( 1491) . and 

|0i) = (1 + 2k(Ax) 2 ) \tJ>) - KX 2 \ip) + 2k{x)x\i(j} - n(x 2 )\?p), (52) 

and its eigenvalue, from ( T4TP and (1511) . is 

Ai = 1 - 2k(Ax) 2 . (53) 

Until now we have made no assumptions about the state (i.e. would-be wave function) 
|^), but in order to have a simple form for the remaining eigenvectors we shall specialize 
from now on to the cases where ip is either an even function of x (i.e. ip(x) = ip(—x)) 
or an odd function of x (i.e. ip(x) = —ip(—x)). Using this assumption in ( |52l) and (1531) 
we have that 

\<p l ) = (1 + k(x 2 )) \^)-kx 2 \^), 

and 

\ x = 1- 2k(x 2 ). 
Next we note that any eigenvector \<f>) must (to order k) satisfy 

P *\<p) = \ij)(m-Kx 2 \ip)(m 

+ 2Kx\ip)(ip\x\(j)) - fi#)(^|x 2 |0) = \\<f>). (54) 

from which it can be seen that any non-zero eigenvector of p K , must be a linear 
combination of x\ip) and x 2 \ip). Using this observation as a clue, it is easy to 
see that there is another eigenvector, |0 2 ), with a non-zero eigenvalue, and that it has 
the simple form 

|02>=3#>, (55) 

with an eigenvalue equal to 2k(x 2 ), which can be seen from ( 1551) and fj4"5]) . 

Since the sum of the eigenvalues we have found so far is 1 (= tr p = tr p K ) and on 
the assumption that p K , like p is positive, it must be that the orthogonal complement to 
the (two-dimensional) subspace spanned by the eigenvectors we have found so far must 
consist entirely of eigenvectors with eigenvalue zero (to first order in k). Thus, to first 
order in k, we have found all eigenvectors and eigenvalues of p K and (to this order) these 
are obviously the same as for p. 

We remark that an alternative derivation of the same result, and confirmation 
that, for sufficiently small k, p K is positive, may be had by noting that (still with our 
specialization to ip either even or odd) to first order in k, (T451) may be written 

p = |(1 — KX 

V))((l - KZV)! + 2k|(2^))<(2^)| 

This is in the canonical form, 

P = y^&nl^n) (lj>n\ 
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except that the vectors which appear, while orthogonal, are not normalized. To remedy 
this, we easily have (to first order): 

||(1 - KX 2 i))\ \ = 1 - K(lp\x 2 lp) — \ — (IK 

(where, for convenience, we resume our abbreviation of (ip\x 2 ip) by a) so that (again, 
all equalities are to first order) the normalized vector is 

(1 + <m)(l — nx 2 il)) — 1 + an — kx 2 i\) 

and also 

| HI | =a 1/2 . 

We can now replace the vectors by their normalized forms at the expense of introducing 
prefactors equal to the square of their norms. Thus 

p = (1 — 2<m)|(l + clk — K,x 2 tp))((l + a,K — Kx 2 ip)\ + 2aK\(a~^xip))((a~^xxl))\, 

from which we can read off the eigenvectors and eigenvalues. We see that these are the 
same as those we obtained (by a rather different route) above. 

9.1. Diagonalizing a tensor product state in one dimension 

Consider the tensor product state 

= |^) ® |^ 2 ) (g, ... |^), 

and the density matrix 

p(x 1: .., x N ; x' 1: x' N ) = .., x N )V*(x' 1: a / JV ) e -«(«i+^+-+*»-*i'-^'--^') 2 

= V( Xl , .. : x N )^*(x' 1: ...,x' N ) (l -k[(x 1 + x 2 + ...+x N ) - (xi +x 2 ' + ... +x N ')} 2 ^j , 

to O(k). (From now on, we shall not always bother to mention that we are working to 
order k and we shall not bother to distinguish between p and the counterpart to what 
we called above p K .) In Dirac notation this is 

p = |^)(^|_ K ( a;i + ... +a ; 2 ) 2 |^)(^| 

+2k(x! + ... + x N )\^f){^\(x 1 + ... + x N ) - /c|*)(*|(a;i + ... +x N f. (56) 

From the previous example we are tempted to hope that there will once again be only 
be two eigenstates of our density matrix, and that they will have a similar form. We 
know that there is only one eigenstate of the density matrix p = |\&)(\&|, namely |\&). 
Thus we expect that there will be an eigenstate close to \^/) with eigenvalue slightly less 
than unity. Let us write our first eigenvector as 

|$i) = |*>+«I*l> (57) 
where = 0. Then 

p|0i) = | - k(xi + ... + a; 2 ) 2 |^)(^| + 2«(xi + ... + x 7V )|^)(^|(a; 1 + ... + x N ) 

-«|*)(*|(a;i + ... + x N ) 2 } + k\V ± )) = (1 - «a)|0i) 
= (l-«a)(|*> + «|* ± )). (58) 
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Equating the 0(k) terms in the above equation yields 

- (xi + ... + xtv) 2 |*) + 2(xi + ... + xat)|^)(^|(xi + ... + x N )\V) (59) 

+ ... + x N ) 2 \^) = -a|tf) + 

In order to proceed and to be able to find the remaining eigenvector we assume that the 
wave functions ipj(x) are either even or all odd for all j G 1...N, bearing in mind that 
1^) = ® ••• ® iV'iv)- That is to say that for each j we have either ipj(x) = ipj(—x) or 
ipj(x) = —ipj{—x). Assuming this and acting on both sides of the above equation with 
(\&| gives 

a = 2((x 2 1 )+... + (x 2 N )). (60) 

Substituting (|60|) into (1581) we are able to deduce \^?±), and hence the first eigenstate 
(by substituting \^f±) into (1571) ). 

|0l) = [1 + K ((Xi 2 ) + (X 2 2 ) + ... + (X N 2 ))] \V) -K( Xl +X 2 + ... + X^) 2 |tf ), 

along with its eigenvalue 

Ax = 1 - 2k ((x! 2 ) + (x 2 2 ) + ... + (x N 2 )) . 
It can be shown that there is only one other eigenvector, 

\<h) = + ... + x N )\^), 

with eigenvalue 

A 2 = 2k (( Xl 2 ) + (x 2 2 > + ... + (x N 2 )) , 
which is similar to the previous example. 

9.2. Diagonalizing the density operator in three dimensions 

Let the wave function ^(x) satisfy the following symmetry conditions 

^(x, y, z) = ip(-x, y, z) = 4>(-x, -y, z) = ip(-x, -y, -z) = (61) 

ip(-x, y, -z) = ip(x, -y, z) = ip(x, -y, -z) = ip(x, y, -z), (62) 
and be such that it can be written as a product of three independent functions, 

ip(x,y,z) = a(x)f3(y)^(z). 
Consider the density operator 

p(x, x') = p (x, x')e- K(x - x ' )2 ~ Po (x, x') - k(x - x') 2 p (x, x') 

where po(x, x') = ^(x)^*(x'). We wish to find the eigenvectors and eigenvalues of p to 
O(k), which means that we want to find 0(x) and A that satisfy 

J d 3 x'p(x,x')<f)(x') = A0(x). (63) 

In Dirac notation our density operator, to first order in k, is 

p = \ip)(tp\ - «x 2 |V>)(V>| + 2kx|V)(V'| - 2«#)(</>|x 2 . 
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To find the first eigenvector we use a similar perturbative argument to the previous one 
to deduce that there must be an eigenvector close to with eigenvalue slightly less 
than unity, i.e. 

01 (x) = ^(x) + /#_l(x), 

or in Dirac notation 

|0i) = |V>) + «hM. (64) 
We want to find \ip±) such that |0i) is an eigenvector of p. To do this we first write 

p\<f>i) = - ^x 2 |^| + 2/o#)(V>|x - K#>(^|x 2 ) + «#±)) 

= (l-o«)(|V) + /c|V±». (65) 

From this, and using the symmetry restrictions we imposed on we obtain \ip±), 

hM = -(x 2 + (x 2 )M. (66) 

By acting on both sides of fl65l) with we find that a = 2(x 2 ). Thus our first eigenvalue 
is 

Ai = 1 -2(x) 2 

and the corresponding eigenvector (obtained by substituting (1661) into (I641 ). 

|0 1 ) = (l- K (x 2 ))|^)- K x 2 |V). 
We claim that the following are eigenvectors of p, 

2 (x) = xip(x); 02 (x) = W(x); 3 (x) = *V(x), 
with respective eigenvalues 

A 2 = 2k(x 2 ); A 3 = 2k(?/ 2 ); A 4 = 2n(z 2 ). 
To show that 02 (x) is an eigenvector we simply show that it satisfies (163]) through 

J <i 3 x'p(x, x')0 2 (x') = y dx'dy'dz'a(x)P(y)^(z)a*(x')P*(y')^*(z')x'a*(x')P*(y')'j*(z') 

-k J dx'dy'dz'(x 2 + y 2 + z 2 - 2xx - 2yy' - 2zz + x' 2 + y' 2 + z' 2 ) 
a(aO/%h(*X (y'h* (z')x'a* (x')(3* (y'h* (z'). 

The symmetry restrictions on ip(x, y, z) imply that the first term in the above equation 
vanishes. Noting this and then expanding out the second term gives 

J d 3 x'p(x, x')0 2 (x') = —Ktp(x, y,z)(x 2 + y 2 + z 2 ) 

x J dx'dy'dz'a*(x')[3*(y')j*(z')x'a(x')[3(y')j(z') — Kip(x, y, z) 

x J dx'dy'dz'a^x^^iy'h^z'^x' 2 + y' 2 + z' 2 )x'a(x')(3(y'h(z') 

+2nip(x, y, z) I dx ' dy ' dz a* (x) f3* (y')"f* (z')(xx' + yy' + zz)x a{x')f3{y') r ){z) . 
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Once again the first two terms of the above equation vanish and we are left with 

which shows that <p2 (x) is an eigenvector. It is easy to show that 03 (x) and 04 (x) 
are eigenvectors by a similar argument. This example, with its artificial symmetry 
restrictions on the wave function ^(x), serves as an indication that there are more 
eigenstates of p when investigating higher-dimensional systems. 

10. Notes 

(i) The problems and puzzles which concern quantum physics in general include: 

measurement Those usually collected under the heading 'the measurement 
problem in quantum mechanics'. 

entropy The puzzling contradiction between the traditional understanding of the 
entropy of a general closed system in terms of (subjective) 'coarse-graining' and 
our understanding, since the work of Bekenstein [33] and Hawking [31] in the 
1970s that the entropy of the specific closed system consisting of a black hole 
(either sitting in and radiating into, an asymptotically flat empty universe or in 
equilibrium with radiation in a suitable box) is something objective - namely 
something which is approximately equal to 'one quarter of the area of the event 
horizon'. 

second law The problem of supplying a microscopic explanation for the second 
law of thermodynamics, i.e. a microscopic explanation for why the entropy of 
a closed system always increases. 

The problems and puzzles which relate specifically to black holes include: 

information loss The 'information loss puzzle', which we shall take here to mean 
the puzzle as to how it can be that, on the one hand, during the process 
of stellar collapse and black-hole formation and then black-hole evaporation, 
entropy presumably actually continually increases while, on the other hand, the 
underlying quantum dynamics of the process is presumably unitary. We remark 
that what makes this a puzzle is the traditional (but incorrect according to [1] 
and the present paper) assumption that (in this context) the physical entropy 
of a quantum mechanical system should be identified with the von Neumann 
entropy of its (total) density operator - the puzzle then arising because this 
quantity, being a unitary invariant, must therefore remain constant under a 
unitary time-evolution and, in fact, if the (initial) state is a vector state, be 
zero (at all times). 

coincidence The problem (raised and solved in pQ) of explaining the coincidence 
that the thermodynamic entropy ^thermodynamic (-^0 °f a (assumed for 
simplicity spherically symmetric, neutral) black hole of mass M turns out to 



Quantum gravitationally decohered quantum mechanics 



45 



have (approximately) the same value - namely l/47rM 2 , i.e. 'one quarter of the 
area of the event horizon' - as the gravitational entropy ^Q^Q^g^awking^^ 
of a Gibbs state of matterless gravity at temperature T when T is taken to be 
the Hawking temperature T = l/8irM. 

Here, what we mean by the thermodynamic entropy ^thermodynamic (^0 ^ s 
the value of S implied by the standard thermodynamic relation 

dE = TdS 

when we identify E with M and T with l/8nM and take S to be when 
M — 0. This is easily calculated to be 4nM 2 . On the other hand, what 

we mean by the gravitational entropy ^Q^on^Hawking^) * s ^ e vame °^ & 
calculated using the standard equilibrium-statistical-mechanics formula 

S = \nZ - (3d\nZ/dp 

from the Gibbons-Hawking |31j Euclidean-quantum-gravity partition function 
Z(f3), at temperature T = 1//3, of a Gibbs state of a matterless <dmb> quantum 
gravitational field in a suitable (spherical) box. In fact, when (cf. [3 lj ) this 
partition function is approximated by its zero-loop value exp(— -f c i ass i ca i) with 
-^classical ^ ne classical action, 4ttM 2 , of the Euclideanized Schwarzschild metric 
of mass M = l/8nT (= P/8tt), the value of S obtained in this way is again 
easily seen to be 4ttM 2 . (We note that the matterless gravitational field here 
is a mathematical construct and, from a physical point of view of course, 
counterfactual. What we call the Gibbons- Hawking partition function for this 
is the quantity one would obtain if one were to remove the terms representing 
matter in the equations in [31J - see the next Note <dmt>) 

Other puzzles for which our hypothesis offers natural resolutions include the 
thermal atmosphere puzzle and a puzzle which we raise here which we call 
the neutron-star entropy puzzle. We shall say what these are and give our 
proposed resolutions of them at the end of Note <EJ>. 

(ii) In an exact unified theory at the Planck energy and higher, one would , of course, 
not expect there to be a clear demarcation between 'matter' and 'gravity'. But at 
energies well below the Planck energy, one expects that there will. By 'matter' here, 
we of course mean everything other than gravity. We remark that our hypothesis 
entails that entropy is a quantity which emerges at low energies and will cease to 
have any meaning at energies around and above the Planck energy. 

(iii) (This Note is really a footnote to Note <B>.) To give more details about the 
connection with the work of Gibbons-Hawking, and to explain what the traditional 
assumptions have to say about our coincidence puzzle, we emphasize that our 
statement of our coincidence puzzle refers to the partition function for a Gibbs 
state of a matterless gravitational field (i.e. 'pure gravity') whereas in [31] - see 
also the review and further discussion in [32] - the partition function actually 
considered is that for a Gibbs state of gravity together with matter fields. It turns 
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out that the zero-loop approximation to this is identical to that for matterless 
gravity discussed in the main text. However the one-and-higher loop corrections 
will involve the matter as well as the gravity. In these (and many other, later) 
references the traditional assumption is adopted that it is the Gibbs state of gravity 
together with appropriate matter which represents the physical total state of a 
black hole in equilibrium in a box and therefore they find the total physical entropy 
to be the sum of the same (matterless-gravity) zero-loop contribution (equal to 
one quarter of the area of the event horizon) discussed here in the main text 
with a term arising from the one-and-higher-loop matter and gravity part of their 
partition function which is supposed in those references to represent a correction 
to the entropy due to the thermal atmosphere of the black hole. It has become 
clear since the early work of Gibbons and Hawking that, due to a quantum- field- 
theoretic divergence which is, in its turn, due to the infinite coordinate-distance 
of the horizon from any point outside the horizon when the appropriate (i.e. 
'tortoise') coordinate is used, this one-loop term is actually impossible to calculate 
within the Euclidean quantum gravity framework without imposing an ad hoc cut- 
off near the horizon and its value is therefore in doubt (and, as far as we are 
aware, this doubt has still not been convincingly removed by any other approach 
to quantum gravity). For the traditional assumptions to resolve the coincidence 
puzzle, one would clearly have to assume, as is in fact implied in the Gibbons- 
Hawking references, that the one-loop 'thermal atmosphere' entropy is only a small 
correction to the zero-loop term, so that the total entropy of the traditionally- 
described black hole in thermal equilibrium in a box will be well approximated 
by what we called in Note <B> ^G^bons-Hawking^)" wom d then formally 

resolve the coincidence puzzle within the traditional assumptions, since, as in 
the resolution on our hypothesis, both the quantities ^Qi^oiisHawking^) ^ or 
T = l/8nM) and S'thermodynairiic(-kf) would again be understood to approximate 
the physical entropy of a physical black hole of mass M — in one case in equilibrium 
in a box and in the other case radiating into empty space. But this is surely a 
very unsatisfactory and unconvincing resolution of coincidence since on the one 
hand, the quantity ^thermodynamic (M) derives from the thermodynamics of the 
thermal atmosphere (which, after all is what would escape as Hawking radiation 
were the walls of the box to be removed) while, on the other hand, the entropy of 
the thermal atmosphere has to be asserted to be approximately zero in order to 
make the resolution work! This is of course related to the thermal atmosphere 
puzzle to which, as we explain in Note <EJ>, our hypothesis also offers a natural 
resolution. Indeed, on our hypothesis, the entropy of the matter part of the thermal 
atmosphere is predicted to be approximately equal to ^Q^ons-Hawking^) which, 
up to one-and-higher loop graviton corrections, is 4ttM 2 . The one-and-higher loop 
graviton corrections themselves, on the other hand, are plausibly comparable in 
magnitude to the entropy of just one of the matter fields making up the thermal 
atmosphere, and therefore might be expected to have a value of the same order of 
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magnitude as 4irM 2 /N - where N is some appropriate value for the 'number of 
matter fields in nature' - and hence plausibly small. 

(iv) In general, given a pure (i.e. vector) state with density operator p = on 
a total Hilbert space, H, which arises as a tensor product, 7i m ® 7i g of two other 
Hilbert spaces ('m' and 'g' here could stand for 'matter' and 'gravity' or could be 
given a more general interpretation) then one can define the two partial density 
operators p m on Ti m (the standard partial trace of p over TC g ) and p g on H g (the 
standard partial trace of p over TC m ) and it is well-known and easy to show that the 
two von Neumann entropies S m = — tr(p m lnp m ) and S g = —ti(p g \np g ) are equal. 
In modern terminology, an alternative name for these latter (equal) quantities is 
the m-g entanglement entropy of the state p. 

(v) Our hypothesis may be understood as a specific variant of the environment- induced 
decoherence paradigm [251 126] namely a variant in which the total state of a closed 
system is taken to be a pure state evolving in time according to a unitary time- 
evolution. But it has three crucial new features, (1) that the gravitational field has 
a privileged status as a permanent piece (see Note <dxTIb>) of the total environment, 
(2) that the entropy which results by tracing over it is regarded as an attribute of the 
total closed matter-gravity system rather than only of the matter, and (3) that this 
entropy is identified with the closed system's physical entropy. With these three new 
features, one obtains, as far as we are aware, for the first time, a definite notion of 
decoherence even for closed systems and an objective definition (the matter-gravity 
entanglement entropy) for the physical entropy of a general closed system, thus 
offering a resolution to the entropy puzzle (see Note <0>). Note of course that 
this resolution entails that one rejects both the option of identifying the physical 
entropy of a closed system with the von Neumann entropy of its total density 
operator Ptotal (which would of course be zero if, as our hypothesis assumes, p is 
pure, i.e. a projector onto a single vector) and that one also rejects the traditional 
attempt to define it in terms of some sort of (subjective) coarse-graining of Ptotab) 
(For this to be a resolution of the entropy puzzle, we of course need to argue 
that in the specific case of a closed system consisting of a black hole, the matter- 
gravity entanglement entropy will have the correct value - i.e. (approximately) 
one-quarter of the area of the event horizon. We postpone our argument for this to 
the next paragraph.) Moreover, by adopting our new general definition of entropy, 
one obtains a plausible mechanism for an entropy-increase result for closed systems 
thus offering resolutions to the puzzles second law and information loss as we 
explained in the main text. 

Turning to the application of our general hypothesis to the subject of black hole 
thermodynamics, our hypothesis goes naturally together with a radically different 
description of quantum black holes from the traditional description. In this new 
description, the total state of a (say spherically symmetric, neutral) quantum black 
hole of mass M in equilibrium with radiation in a suitable box is understood to 
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be a pure state of a quantum-gravitational closed system and not, as would be 
traditionally assumed, a Gibbs state at the Hawking temperature. However, this 
pure total state is expected to possess, in the relevant low energy description, 
just the right sort of matter-gravity entanglement for its partial trace p ma tter over 

n 

gravity ^° resemble a Gibbs state of the matter fields at the Hawking temperature 
T = 1/ 8ttM in the presence of a background black hole of mass M, and also for 
its partial trace, Pg rav ity, over 7i ma tter; t° resemble a Gibbs state of (matterless) 
gravity at the same Hawking temperature. In fact, it is reasonable to assume, 
the latter state may be identified with the Gibbs state of a matterless quantum 
gravitational field - whose partition function may be calculated by the Gibbons- 
Hawking Euclidean path-integral method summarized in Note <0> by removing the 
matter terms in [31 J. (See Note <dml>.) We remark that of course this extension 
of our hypothesis entails that the appropriate 'theory of everything' will turn out 
to have states which, at suitably low energies admit of such a description with 
appropriate matter multiplets etc. 

We are now in a position to complete our resolution of the entropy puzzle, for, by 
the result mentioned in Note <JTvb> . the matter-gravity entanglement entropy of our 
black hole state must equal the von Neumann entropy of Pg rav ity and therefore, 
in view of the second assumed resemblance mentioned above, must be close to 
'one quarter of the area of the event horizon' (47rM 2 ). Moreover we see that, as 
a byproduct of this discussion, our hypothesis has led to a natural resolution to 
the coincidence puzzle since, having identified the matter-gravity entanglement 
entropy with the physical entropy of our black hole, we would also expect it to be 
the correct '5" in the thermodynamic relation l dE = TdS\ 

Our hypothesis also offers a natural resolution to some other puzzling features of the 
traditional understanding of black hole entropy. Thus, previously, it was unclear to 
what extent the entropy of a black hole (say in equilibrium with radiation in a box) 
resides in the quantum geometry of the hole itself, i.e. in the quantum gravitational 
field, and to what extent it resides in the matter-fields, i.e. in the thermal 
atmosphere. This is essentially the thermal atmosphere puzzle discussed e.g. 
in [5]. It is particularly puzzling because, on the traditional understanding, it 
was possible to make arguments (cf. Note <dmt>) to the effect that the entropy 
all resides in the gravitational field but one could also (cf. also e.g. the recent 
review by Don Page (35]) make counter- arguments to the effect that it all resides 
in the thermal atmosphere. As a further byproduct to the discussion above, we 
see that our hypothesis also offers a very neat solution to this puzzle. For if we 
identify 'the entropy of the gravitational field' with the von Neumann entropy of 
^gravity an d 'the entropy of the (matter part of) the thermal atmosphere' with 
Pmatter an d the entropy of the total state with the matter-gravity entanglement 
entropy, then, by our resolution to entropy (combined with the result mentioned in 
<Hvl>) we see that our hypothesis entails that all these three quantities are actually 
identical. In particular it turns out on our hypothesis that both statements 'the 
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entropy resides in the gravitational field' and 'the entropy resides in the thermal 
atmosphere' (which traditionally were regarded as mutually exclusive possibilities) 
become simultaneously correct statements! 

Another puzzle in the traditional understanding (which we point out here) was 
that while a black hole (formed say from neutron star collapse) is traditionally 
understood to have a very large non-zero entropy, the entropy of a neighbouring (not 
necessarily stationary) state of quantum gravity which consists just of a neutron 
star (on the verge of collapsing) would, traditionally, be assigned a zero entropy. We 
call this the neutron star entropy puzzle. With our definition (matter-gravity 
entanglement entropy) any state of any closed system would be expected to have a 
non-zero entropy and the entropy of such a neighbouring neutron-star state could, 
quite plausibly, have a value which neighbours (but is a bit less than) the entropy 
of the black hole it is on the verge of collapsing to. (For preliminary estimates of 
the entropies, on our hypothesis, of certain simple model matter-gravity systems 
in the Newtonian limit, see [2j. These have very small values, but, as discussed in 
Section [7J the Newtonian results are expected only to be "the small non-relativistic 
tip of a large relativistic iceberg"!) 

What our hypothesis tells us about the puzzle measurement is of course the 
subject matter of the rest of this paper. 

Finally, we remark that one might regard the fact that our hypothesis allows a more 
satisfactory definition for the entropy of an open system than traditional ideas do 
as further evidence for the correctness of our hypothesis. See Note <lxii[>. 

(vi) As explained in [2], e ~ D(ai ' a2) is the value of the inner product (gi\g2), in a 
suitable Hilbert space, Hg rav ^ y , for the linearized quantized gravitational field, 
between two (non-radiative) state vectors, gi,g2 £ "^gravity which are the quantum 
representations of the static Newtonian gravitational field due to our ball at rest 
with its centre of mass located at ai, a 2 respectively. (As discussed in [2], (<7i|<7 2 ) 
turns out to be real.) To understand how the formula (j3J) arises, it is helpful to first 
consider a would-be wave function, ip, which consists of two sharp peaks localized 
around, say ai and a 2 ; we shall write 

ip = ci^i + c 2 ip2 (67) 

where c\, c 2 6 C, ip\ is a normalized wave function consisting of single sharp peak 
centred on ai and ip2 is a normalized wave function consisting of a single sharp 
peak centred on a 2 . 

In this case, it is envisaged that a full quantum gravitational description is given 
by a vector in 7i ma tter ® ^gravity which takes the form 

* = cxipt <g> gi + c 2 ip2 ® 92- 

p, obtained by tracing |^ r )(^ r | over Hg rav ^y, is then clearly given by 

p = IcilViX^il + cic* 2 (g 2 \gi)\ipi) {4>2\ + c*c 2 (s<i|s' 2 )|V ) 2)(V'i| + M 2 !^)^! (68) 
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and a little thought should suffice to see that, (a) the sharper the two peaks in the 
would-be wave function, the closer this will be, in its position space representation 
to (J3]), (b) for a general would-be wavefunction, one may argue that the formula (jSJ) 
must still hold by regarding it as a limiting case of a wave function with multiple 
peaks. 

(vii) The generalization of the Gaussian asymptotic regime to N balls would be 
relevant e.g. to a superposition of two configurations, each involving N balls, 
such that, while the balls within each configuration are possibly well separated, 
the configurations themselves differ only slightly so that each ball in the second 
configuration is displaced by much less than its radius relative to its position in the 
first configuration. We remark that the formula easily generalizes to the case of N 
balls of equal radius R but possibly differing masses in which case D takes the form 

9 ( M tW^)(^m-4m) 2 

where M^gj is sum of the N masses and x C m and x' cm the centres of mass of the 
unprimed and primed configurations which constitute the arguments of D. We also 
remark that the generalization to the many-ball case of the logarithmic asymptotic 
regime has 

N N / 1 / II / | \ -\2MiMj 

exp(- J D(x 1 , . . . , x w ; x' 1; . . . , x^)) = \\ ] [ , , i' ~ */ (69) 

where it is to be understood that, in the cases I = J, the terms in the denominator 
|xj — xj| |xj — x'jl are to be replaced by R 2 j (This is an easily obtained generalization 
of the formula given ('for simplicity') in [2] in the case where all the Mi and all the 
Rj are the same. The robustness of the latter formula under changes in shape and 
graininess of the bodies is discussed in [3]. 

(viii) The rule (|HJ) is an automatic consequence (see Note <JvTl>) of taking the total state 
for a would-be wave function ip('x), say with a single sharp peak at x = a, to be 
ij; <S> g where g is the (non-radiative) quantum counterpart to the static Newtonian 
gravitational field of our ball when its centre of mass is located at a - irrespective 
of whether the would-be wave function itself is static or not. Obviously, it amounts 
to neglecting the radiative part of the gravitational field. This approximation and 
its generalizations to the many-ball case etc. are however expected to be excellent 
as long as the classical gravitational radiation from the relevant systems would 
be negligible. See Note < bcxvib > for a quantitative estimate on how good this 
approximation is likely to be. 

(ix) Ignoring some technicalities, the GKS/L form (often known as 'Lindblad form') 
for a master equation holds if and only if the master equation integrates up, for 
positive times, to a 'semi-group of completely positive maps' (sometimes known 
as a 'quantum dynamical semigroup') on the relevant space of density operators 
- see [6], [7J. ([6] treated such semigroups of matrices and [7J of bounded 
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operators.) A prototype (see Section E]) is the BLP master equation (1551) [27\ 
p = —i[H,p] — c[x, [x, p]), c > 0, for a one-particle density operator. 
If one has a total Hilbert space, H, which arises as a tensor product of form TC m ®7~tg, 
and if the dynamics of a total state vector ty(t) G Ti evolves according to a (unitary) 
Schrodinger time-evolution for some Hamiltonian, H, then, if the time-zero state 
vector ^(O) takes the product form <p m ® ip g , then one knows (|6|) that for any 
(positive or negative) time t, the partial trace, p m {t), of p{t) = |\E'(t))(\I/(t)| over 7i g 
will arise by the action on the (pure) time-zero density operator p m (0) = |</> m )(</> m | 
of a 'completely positive map' (which will of course depend on ip g ). (See e.g. also 
[6] for the definition of 'completely positive'.) 

Moreover, for positive times and if TC m represents now, say, 'charged stuff' and H g 
now, say, 'radiation', and if H is such that the 'produced radiation' tends to be 
emitted from, and fly away from, the 'charged stuff', then one expects, in line with 
general expectations (see [6] and references therein) that the time-evolution will be 
approximable by the action on p m (0) of a semi-group of completely positive maps 
whose exact solutions obey a master equation of GKS/L form. (This semigroup 
will extend to act on arbitrary - i.e. not-necessarily-pure p m (0)) We remark that, 
although this is not always mentioned or emphasized, one of course equally expects 
the time evolution for negative times to be similarly governed by a master equation 
of 'anti-GKS/L' form (i.e. a generator of a completely positive semigroup in the 
variable t' = —t) and this will obviously be the case if the Hamiltonian H and initial 
state \1/ (0) mentioned above are both time-symmetric in obvious suitable senses. 
Benatti and Narnhofer [36j have given necessary and sufficient conditions for the 
von Neumann entropy of a density operator evolving in time (for positive times) 
according to a master equation of GKS/L form to increase monotonically. One can 
easily check that their conditions hold for BLP and we mention in passing that it is 
shown in [36J that they hold for the time-evolution of GRW |13j . (The same result 
with S replaced by S\ - defined in Section [6]- is demonstrated in Section [6] for BLP 
and in [13] itself for GRW.) When a total unitary time evolution as in the previous 
paragraph leads to an approximate GKS/L-form master equation governing the 
behaviour of p m at positive times which satisfies Benatti and Narnhofer's conditions, 
then one clearly expects entropy to increase for positive times. Although this is 
again not always mentioned or emphasized, one also expects (and again of course 
this must happen if the total dynamics and time-zero state has a time-reversal 
symmetry) the entropy of p m will increase for increasingly negative times. We shall 
call such an overall time-behaviour for the von Neumann entropy of p m as t ranges 
over the whole real line a 'two-sided entropy increase' result. 

(The above two paragraphs are particularly relevant to the discussion in Section 

El) 

(x) To see that the solutions to (GO) and ( TTOi) specified by the formulae (IE]) and (J6j) 
consist of density operators at all times, we note that it is clear from the fact that 
-D(O) = 1 that p will have unit trace. It is also clear from the way that ([3]) and 
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(|SD etc are derived that p(t) will be a positive operator at all t, and this can be 
shown to hold also when D in (jSj) is replaced by either its Gaussian or logarithmic 
asymptotic forms. 

The failure of the master equations (Q and (flOl) to have GKS/L form can be traced 
back to the fact that the modes which are traced over in obtaining Pphysical ( see 
Note <JvTb>) are, in the Newtonian case, non-radiative and, in contrast to what was 
envisaged for 'radiation' in our discussion of the approximate validity (for positive 
times) of master equations of GKS/L form in the situation envisaged in Note <drxb>. 
instead of 'flying away' from our matter, they are 'dragged around with' (or one 
might say 'slaves to') the matter. (Cf. the last paragraph of our discussion of the 
Penrose experiment in Subsect ion 14 . 1 1 and Note < ixxvi > and also (especially the last 
paragraph in) Section [6j 

Below, we shall call time-evolutions on spaces of density operators of form described 
in the paragraph containing equation (|Sj) 'multiplicative decoherence models' 
(MDM) and master equations such as (E]) and (TTOT) 'MDM master equations'. 
It is interesting to compare and contrast MDM master equations with master 
equations of GKS/L form (see Note <drxb> and Section [6]). MDM master equations 
have solutions for all times, t, ranging over the real line while GKS/L master 
equations only have solutions for positive t. On the other hand, for positive 
t, GKS/L master equations have solutions for arbitrary initial density operators 
while MDM master equations only have solutions for very special initial data (but 
including of course initial data of form (jHJ) for equation (jHJ) and (jSJ) for equation 
(TlO|) etc.) What goes wrong for GKS/L equations for negative times, and for MDM 
master equations for wrong initial data is that one can find solutions in the class 
of unit-trace trace-class operators, but they will fail to stay in the class of density 
operators because they will fail to satisfy positivity. Related to this, one can prove 
(BSK unpublished), for our Gaussian model master equation ( flQj) . the theorem: 

The time evolute after an arbitrarily small positive time 5t of a pure density operator 
p = where ip is even, fails to be positive except when ip takes the Gaussian 

form ip = e~ cx2 with Rec > and Imc < 0. (In the special case with c = Ak, 
p = evolves into itself for all later times.) 

If, however, we compare and contrast MDM models, not with time-evolutions 
generated by master equations of GKS/L form, but rather with time-evolutions 
such as the sort described in <drxb> which, as we discussed in that note, are expected 
to be approximately of GKS/L form for positive times and of anti-GKS/L form for 
negative times then, as far e.g. as the behaviour of entropy is concerned, the 
two sorts of models can have qualitatively similar behaviour. Namely, in both 
cases (depending on the details of the dynamics) one can have a two-sided entropy 
increase result in the sense explained in <drxb>. We show such a two-sided entropy- 
increase result for our Gaussian (MDM) model with a free Hamiltonian (and to 
first order in k) in Section [6j 
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(xi) As we shall discuss further in Section [U we expect that our conclusions will also 
apply if the probe particles are relativistic (e.g. photons). 

(xii) The purpose of this Note is to sketch a natural extension of the hypothesis of [1] 
(see subsection II. II of the Introduction and also Notes <0> and <EJ>) to deal with 
open systems. In this extension, the specification of an open subsystem of some 
given closed system is taken to correspond to a particular way of expressing the 
Hilbert space 7i ma tter as a tensor product 

^matter = ^matter, system ® ^matter, environment 
so that, in view of (CD), the Hilbert space, T^total f° r the total system will arise as 
a triple tensor product 

^total = ^matter, system ® ^matter, environment ® ^gravity 
For a given Ptotal (which, see ([2]), will, according to our hypothesis, take the form 
|^ r )(^ r |) of the relevant total closed system, the density operator, p ma tter, system 
describing the partial state of the matter belonging to the open subsystem is then 
declared to be the partial trace of p total over ft matt er, environment ® ^gravity and 
the physical entropy of the open subsystem is declared to be the von Neumann 
entropy of Pmatter, system- All this may be represented by the schematic rectangle- 
picture in Figure 4a where one may regard the vertical dividing line as separating the 
matter degrees of freedom which belong to the 'open subsystem' (marked 'matsys' 
on the figure) from those which belong to its 'environment' (marked 'matenv' on 
the figure) while the horizontal line separates all matter from gravity (indicated by 
the region marked 'grav'). In the sense that it is always part of what is traced over, 
one might say that gravity is regarded as a permanent part of the environment so 
that the horizontal line is fixed, while the vertical line is slideable. Sliding it to the 
right corresponds to gradual enlargement of our open subsystem of interest until, 
in the limit as it coincides with the right hand boundary of the lower rectangle, the 
open system will approach the total closed system and its entropy will approach the 
entropy of the total closed system which, of course, will still be non-zero because all 
of H 

gravity continues still to be traced over in the definition of that. Indeed, if our 
total closed system contains a (say stellar or galactic-centre sized) black hole, then 
this limiting value of the entropy will be very large. As we indicate in the schematic 
graph Figure 5a, it seems reasonable to expect that the entropy of ever larger (say 
nested) subsystems will increase monotonically with the size of the subsystems. 
All this is in contrast to, and we feel, more satisfactory than, the situation if one 
were to attempt a definition for the entropy of an open subsystem of a given closed 
system on traditional ideas (i.e. on the assumption that the gravitational field 
due to the matter of the subsystem is regarded as belonging to the subsystem) by 
defining it to be the von Neumann entropy of the partial trace of the total density 
operator over the Hilbert space for the total environment of the subsystem, (and 
assuming the total density operator to be pure). In other words, by defining it 
to be the entanglement entropy of the subsystem with its total environment (and 
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Fig. 4a 




Fig. 4b 



Figure 4. Schematic diagrams contrasting our approach to open systems (Fig. 4a) 
with that on the traditional 'environment-induced decoherence' paradigm (Fig. 4b) 



assuming the total density operator to be pure). This might be represented by the 
schematic rectangle figure 4b. The entropy, thus defined, for a given open subsystem 
would then necessarily (cf. Note <Hvt>) equal the entropy of its environment and 
therefore, as one considers ever larger nested subsystems which eventually approach 
the total closed system, as indicated in figure 5b by sliding the vertical dividing line 
to the right, the adage 'what goes up must come down' would necessarily apply and 
the entropy would approach the value zero! So a typical graph of entropy against 
size of region would, in contrast to Figure 5a look like that sketched in Figure 5b. 




size of region size of region 

Fig. 5a Fig. 5b 



Figure 5. Schematic behaviour of entropy against 'size of open system' for our 
approach to open systems (Fig. 5a) contrasted with the corresponding behaviour 
on traditional ideas (Fig. 5b) 

Finally, we remark that, in defining the density operator of any open system, the 
result of tracing over W ma ^er, environment ® ^gravity * s °f course equivalent to first 
tracing over 7ig rav ity an< ^ then tracing over ?i ma tter, environment- So our proposal 
is equivalent to a modification of the above-described approach to defining entropy 
on traditional lines - i.e. as pictured in Figure 4b - in which one identifies 'system' 
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as 'matter system' and identifies 'environment' as 'matter environment' but, instead 
of postulating a total state which is pure, one postulates a total state on the full 
matter system which is equal to the (mixed) trace of our matter-gravity Ptotal over 
^gravity- 

(xiii) In this note, we shall add a number of further remarks about our resolution to the 
Schrodinger Cat puzzle discussed in Sections [T] and 

We begin by remarking that, in the context of Schrodinger Cat-like situations, 
something similar to our mechanism whereby a pure state such as (1111) gets changed 
to a mixed density operator like fl68l) or (approximately - see below) like ffT2l 
occurs, of course, in the general context of 'environmental induced decoherence' 
(see e.g. [25] and [26] and references therein). However there are a number of 
important differences. First, the 'environmental induced decoherence' mechanism 
is only capable of bringing about the replacement of a pure by a mixed state in 
the case of an open system, whereas in the theory of pQ and [2], by having a 
preferred environment (i.e. gravity) and declaring that this must always be traced 
over in principle, such a replacement of pure by mixed states occurs for closed 
systems. Secondly, something similar to our proposed general interpretation (see 
around the passage labelled (B) in Section CD and the beginning of Section [5]) of 
our physical density operators in terms of events which happen is usually implicity 
or explicitly assumed when discussing the interpretation of density operators of 
form (fl2l - but it usually seems to be assumed that the events which can happen 
should be identified with one-dimensional diagonalising subspaces of the relevant 
density operator (and that the probability with which they happen is the associated 
eigenvalue). In the case of ffl2l this will not change the interpretation as long as 
| ci | 2 7^ But in that special case, on the latter usual interpretation, an ambiguity 
occurs since the eigenvalue ~ is degenerate. The existence of this ambiguity in this 
usual interpretation is sometimes taken as a reason to question whether a formalism 
based on density operators can really solve the conceptual problems of quantum 
mechanics. However our interpretation (B) sidesteps this criticism by entailing 
that, in this same special case, there is only one possible event which can happen 
(with probability 1), namely that corresponding to the subspace spanned by ipi 
and ip2- This is free from ambiguity although it is admittedly somewhat strange. 
However, the strangeness is perhaps mitigated in that, in a typical preparation of 
a would-be Schrodinger Cat-like state, such as ( |TTi) . the quantities c\ and c 2 would 
be expected to vary in time, and the equalities |ci| 2 = ~ = |c 2 | 2 to only hold for 
a single moment of time, and hence our interpretation would entail that there are 
two possible events, each evolving in time, then fusing into one possible event for 
a fleeting moment before splitting into two (evolving) possibilities again. The ball 
(or the cat!) is still predicted to be one thing or another except at that moment 
and at no time is the statement about the possible events which can happen (and 
the probabilities with which they will happen) ambiguous. 
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Finally we remark that (gi\g 2 ) i n (1681) can be small but never actually zero and 
therefore the density operator (|12p can only be an approximation to fl68l) in Note 
<dvil>. In fact, in general the subspaces belonging to our two events will clearly be 
one-dimensional subspaces spanned by (normalizable) vectors of the form 

0i = kipi + s*ip2 and <fi 2 = -sipi + k*ip 2 (70) 

with |fc| 2 + |s| 2 = 1 and p will take the form 

+ J B|0 2 )(0 2 | (71) 

where A and B are positive numbers with sum 1 and the values of A, B, k and s 
will be determined (by equating ( 17T1) with (l68il ) in terms of |ci| 2 , \c 2 \ 2 and {gi\g 2 ). 
When, the ball is of mass much smaller than the Planck mass and/or the two ball 
states are close together, (gi\g 2 ) will be close to 1, A will be close to 1, B close to 
zero, k close to C\ and s close to c 2 so there will be one 'event' with a probability 
close to 1 which is close to the subspace spanned by ip of ( fTTT) while, at the other 
extreme, returning to the case of almost complete decoherence, as our ball gets 
larger and/or as the two ball-states get further apart, k will be close to 1 and s 
close to zero, and A will be close to |ci| 2 and B will be close to \c 2 \ 2 and k will 
get closer to 1, and s closer to zero. We feel that our statement that this "is just 
what one would hope for from a resolution to the Schrodinger Cat puzzle" remains 
equally arguable with this exact description of the event-subspaces. 

(xiv) As far as we are aware, an interpretation in terms of 'events which happen' along 
the lines of that proposed here in (B) has not been considered previously in the 
context of models for modified quantum mechanics which are couched in terms 
of time-evolving density operators. However, something similar to (but not quite 
identical to) our proposal seems often to be implicitly assumed when considering 
the physical interpretation of density operators which arise in contexts such as 
proposed resolutions to the Schrodinger Cat puzzle - see Note < lxiiil >. 

(xv) Here, we reserve the term 'collapse models' for theories such as the 'GRW model 
[T3] when these are formulated in terms of time-evolving density operators. For a 
given such collapse model, (and provided, as is e.g. the case for the GRW model, the 
time-evolving density operator arises from a semigroup of completely positive maps 
- see Note <dixb> - acting on some initial density operator) it is now understood that 
there are typically several different ways of representing the dynamics in terms of a 
time-evolving stochastic wave function - each constituting what is (since the term 
was coined around 1993 by Carmichael [15]) now called a different 'unravelling' 
(see e.g. [II]) of the model. An early example is the 'usual' unravelling of the 
GRW model which actually predates the use of the term 'unravelling'. This 
was formulated in a simple explicit way by Bell in [16]. For a given collapse 
model, a choice of unravelling often suggests another way to that proposed here 
of interpreting the collapse model in terms of 'events' which 'happen'. For the 
'Bell' unravelling of the GRW model, this is explained in [16]. The proposal made 
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here for how a notion of 'events which happen' may be obtained from a given time- 
evolving density operator is offered as an alternative to such interpretations (see 
Note <hrivl >). It involves a more abstract notion of 'events' to the 'collapse centres' 
of Bell/GRW. It may suffer from a number of problems, but it is free from the 
problem of the ambiguity of unravellings. In any case, the notion of unravellings is 
only available for time-evolutions which exactly arise from semigroups of completely 
positive maps - i.e., cf. <dixl>. which obey master equations of exact GKS/L form. 
Hence, unravellings are anyway unavailable to us since the master equations, ( 1421) 
and (H3j) . of relevance to us are not of this form. (At best, (|43|) is approximately of 
this form (for positive times). For similar reasons, unravellings are unavailable for 
more traditional situations where one has a time-evolving density operator arising 
from some sort of environment-induced decoherence ([25], [26]) since such time- 
evolutions also do not constitute semigroups of completely positive maps (these 
again only arising, if at all, after some sort of mathematical limiting procedure). 
On the other hand, the events interpretation we propose here would be available. 

(xvi) It is of course a tricky philosophical issue just what one means by 'probability' here. 
(Some common views about what probabilities are make reference to observers!) 
We will not attempt to say anthing original about this issue here, but we note that 
more or less the same issue is present with notions of events such as that proposed 
by Bell for the GRW model (see Note <|xyj>). A recent article which addresses the 
issue in this latter context is [T7] . 

(xvii) The statements about the way the pencil gets reflected off the mirror can easily be 
verified with a simple conservation of momentum argument making the usual small- 
angle approximations. In explanation of the various assumptions: (a) the diameter 
of the pencil is assumed to be much bigger than S in order to ensure that the angle 
of widening of the pencil due to spreading of the wave packet will be much less than 
the angles between the various beam components after the reflection. (There will 
also be an initial broadening of all the reflected pencils by an amount a/2<5 because 
of the uncertainty in the location of the mirror.); (b) the diameter of the mirror 
has obviously to be large enough for the pencil to be entirely reflected by it, no 
matter where along the interval (0, 25) of the x-axis the centre of mass of the bead 
is located; (c) the mass of the probe particles is taken to be much smaller than M 
so that reflection by the mirror will consist simply in reversal of the perpendicular 
component of momentum in the centre of mass frame; (d) the assumption P ^> p 
was made simply so that we can use the usual small angle approximations to obtain 
the uniform spacing in the angles of the reflected beam components. 

(xviii) A similar declaration to that considered in the discussion motivating our corrected 
pragmatic interpretation in Section H] has been made before in the context of other 
approaches to the measurement problem in quantum mechanics. Cf. e.g. the 
discussion around the following passage in the article [TT] by John Bell concerning 
what might constitute what Bell calls a 'beable': "Not all 'observables' can be given 
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beable status .... What is essential is to be able to define the positions of things, 
including the positions of instrument pointers or (the modern equivalent) of ink on 
computer output." 

(xix) One could e.g. presumably undo the measurement of the bead momentum inserting 
a suitable lens-like device between the probe particle and the screen so that, however 
its wave function reflected off the mirror, only one spot forms on the screen. 

(xx) The version described here is similar to the 'FELIX' experiment described in [21] and 
to the 'suggested space-based experiment' described in Appendix 2 of [22]. In the realistic 
versions of the experiment described in these references, Penrose envisages the photon 
to be an X-ray photon, and the movable mirror to be a Mossbauer crystal. Penrose 
and collaborators have also proposed a ground-based experiment (see again Appendix 2 
in [22] and [9]) in which the photon and movable mirror are optical and the photon is 
stored in suitable optical cavities. The details of the ground-based experiment (see [9]) 
are slightly different but it is easy to see that our discussion here and in the following 
subsection could easily be adapted to this latter experiment and the conclusions we draw 
here would survive in essentially unaltered form. We remark that what is difficult about 
these experiments (and is discussed in the references just given) is of course to eliminate 
all causes of decoherence other than the possible cause that the experiments are designed 
to test. 

(xxi) We have made some extra implicit assumptions and simplifications here. In particular, 
the state of the photon/movable- mirror-cum-oscillator system at time t m might well be 
more accurately modelled by an entangled state, rather than the single tensor product we 
have assumed here. Without extending our analysis to include such a possibility, we note 
that we are still assured by our position measurement theorem that, irrespective of how 
we model it, as long as the Penrose experiment ends in a type of position measurement 
(and assuming that photons may be treated for this purpose as non-relativistic particles) 
the prediction of the theory of [2] will still be identical with the prediction of standard 
quantum mechanics. 

(xxii) We should really have allowed for a relative phase in the two equations (|28p . (|29p . 
However, for simplicity we assume the phase happens to be 1 and that, with this phase, 
one predicts zero detection rate at the detector. If either of these were not the case, one 
could e.g. imagine modifying the experiment by inserting a suitable phase-shifting device 
in one arm of the interferometer in order to continue to predict zero detection rate at the 
detector. 

(xxiii) In his 1996 article, [20J, Penrose explains that what he actually envisages is a radical new 
theory of 'quantum state vector reduction' and he emphasises that the technical form that 
such a theory will take is not yet known. Rather, some qualitative and semi-quantitative 
arguments are given and the direction in which such a theory may lie, indicated. In our 
discussion of the Penrose experiment in Subsection 14.11 we perhaps should have written 
"Now what is envisaged by Penrose, to the extent that it may be represented within 
current quantum mechanical formalism involving density operators etc. . . .". 

(xxiv) Interestingly, we easily see that the partial state of the photon at time t m is predicted to 
be the same as that of (|30p by all three theories: 




(However, in the Penrose experiment, no direct measurement is made of this state.) 
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(xxv) We remark that it is of course the equality of g(tf) with g(to) (or more generally the 
fact that, at any time, t, g{t) depends only on the matter-configuration at time t) in the 
non-relativistic case which is responsible (cf. Note < jviiit >) for the form of (the relevant 
generalization to the photon/movable- mirror-cum-oscillator system of) equation ([8]) which 
(together with the relevant version of equation (|27p ) is responsible for the validity of our 
Position Measurement Theorem (when applied to the Penrose experiment). 

Further insight into the reasons for the recoherence may be had by considering the simple 
model closed system consisting just of a movable- mirror-cum-oscillator as in the Penrose 
experiment, but uncoupled from any photons etc. The centre-of-mass motion of such an 
oscillator is of course modelled in standard quantum mechanics as a standard quantum- 
mechanical harmonic oscillator. Any (Schrodinger picture) wave-function solving this 
will, of course, be periodic in time and taking the would-be density operator, po(t), in 
(a one-dimensional version of) ([8]) to be \tfj(t)){^(t)\ for such a periodic would-be wave 
function, it is clear, from ([8]), that p will also be periodic and hence alternately decohere 
and recohere. This is in contrast with the typical prediction of a 'collapse model' - see 
Section [U] and Notes <dixb> and <[xp> for further discussion. 

(xxvi) By Einstein's quadrupole formula, we would expect the energy radiated per unit time 
to be given by a number of order 1 times m^i^cu® (if we restore c and G, this should be 
multiplied by G/c 5 ) where m is a relevant mass (assuming the mass of the 'oscillator' 
part of the movable-mirror-cum-oscillator [the 'cantilever' in [9]] can be neglected, then 
this would be the mass of the movable mirror and its mounting) and t its amplitude, 
and to its frequency of oscillation. Thus the expected number, -Wgraviton' of gravitons 
emitted per cycle would be a number of order 1 times m 2 l ui (times G/(hc 5 )) (times 
2tt). Putting in the values (similar to or larger than those in the experiment suggested in 
|9j) m = 1CT 8 gram, I = 10 -11 cm, to = 3 x 10 3 sec -1 , this gives an -^g rav iton °^ aroun d 
10 -80 which is utterly negligible! 

We remark that there will actually inevitably be another source of decoherence due to a 
form of radiation by the oscillating mirror, namely the radiation by the movable mirror 
of photons at around the oscillator frequency due to the effect of the oscillating boundary 
condition imposed on the quantized electromagnetic field by the oscillating mirror. Ford 
and Vilenkin [23] have calculated that, for such an oscillating mirror, the energy per unit 
time radiated by this mechanism for a scalar analogue of the photon is given by (note that 
a missing factor of 2tt in [23] has been inserted here) E = 72 Q n '2 ( 2 u>®A (times fo/c 4 ) where 
A is the surface area of the mirror, and we shall assume that the same formula, multiplied 
by 2 because of the photon's two polarization states, holds true for true photons. (It is 
not actually known whether or not this assumption is correct, but it will surely give the 
correct order of magnitude.) The expected number, -^photon' of photons emitted per 
cycle will then be j^l^u^A (times 1/c 4 ). With the same values assumed above for £ 
and to and with a value for A (again taken from the suggested experiment [9j ) of 10~ 6 
cm 2 , this gives an -^p^^^ of around 10~ 39 which is vastly bigger than -^g rav iton 
of course, for the purposes of the Penrose experiment (in the version proposed in [S]), 
reassuringly, still vastly smaller than anything which would need to be taken into account 
in a practical analysis of the experiment. 

(xxvii) The discussion of the relationship between the notion of parity in the naive pragmatic 
interpretation and the notion of parity in our events interpretation and, in particular, the 
coincidence explained after equation (|36p may clearly be generalized as we now explain: 
It seems natural to introduce the following definitions and terminology: For any p to 
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which one intends to apply our events interpretation, we say that a self adjoint operator, 
B, is a beable for p if it commutes with p and if all spectral subspaces of p belonging to 
non-zero eigenvalues of p are also spectral subspaces of B. In other words, a beable for p 
is a self-adjoint operator which commutes with p and, on the orthogonal complement of 
the kernel of p, is a function of p. We then say that an eigenvalue of B on a given event 
(i.e. spectral subspace) of p is called the value of the beable B when that event occurs. 

With this definition, it is obvious that, for a given p, a similar coincidence to that remarked 
upon after equation ([36]) will occur quite generally between the expected value of a beable, 
B, in the sense of the sum over all possible events of the probability of each event which can 
occur multiplied by the value of the beable B when that event occurs and the expectation 
value of B in the conventional sense of tr(pB). 

We have taken the word 'beable' from the writings of John Bell - see e.g. where Bell 
writes: 

"The beables of the theory are those elements which might correspond to elements of 
reality, to things which exist. Their existence does not depend on 'observation'. Indeed 
observation and observers must be made out of beables. v 

However, we should caution that the candidate for a beable being proposed here differs 
from various candidates that Bell himself proposed (see e.g. Note < jxviiil >) and, in 
particular, we would remark that, in contrast to other proposals for beables, in our 
notion, what is a beable at a given time depends on the state of the system at that time. 

A couple of interesting features of our notion of beable are worth pointing out. First if, in 
the case of an n-dimensional Hilbert space, p is 1/n times the identity operator then the 
only beables will be multiples of the identity operator. This applies, for example to the 
Schrodinger Cat state (fT2j) when |ci| 2 = | = | C2 1 2 - (See also the discussion of this special 
situation in Footnote < jxiiil >.) At the other extreme, if p arises as the projector 
onto some vector ip, then any self-adjoint operator which commutes with p is a beable. 
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